To appear in Ap. J. Suppl. 



An Empirical Ultraviolet Template for Iron Emission in Quasars as Derived 

from I Zw l 1 

M. Vestergaard 2 ' 3 ' 4 and B.J. Wilkes 2 
ABSTRACT 

We present an empirical template spectrum suitable for fitting and subtract- 
ing/studying the Fen and Fein emission lines in the restframe ultraviolet spectra of 
quasars and active galatic nuclei, the first empirical ultraviolet iron template to cover 
the full range of A1250-3090A. Iron emission is often a severe contaminant in optical- 
ultraviolet spectra of active galactic nuclei and quasars. Its presence complicates and 
limits the accuracy of measurements of both strong and weak emission lines and the con- 
tinuum emission, affecting studies of line and continuum interrelations, the ionization 
structure, and elemental abundances in active galaxies and quasars. Despite the wealth 
of work on modeling the quasar Fen emission and the need to account for this emis- 
sion in observed quasar spectra, there is no ultraviolet template electronically available 
to aid this process. The iron template we present is based on Hubble Space Telescope 
spectra of the Narrow Line Seyfert 1 galaxy, I Zwicky 1 (I Zw 1, z = 0.061). The intrinsic 
narrow lines (> 900 kms -1 ) of this source and its rich iron spectrum make the template 
particularly suitable for use with most active galactic nuclei and quasar spectra. 

The iron emission spectrum, the line identifications, and the measurements of ab- 
sorption and emission lines are presented and compared with the work of Laor et al. 
Comments on each individual line feature and the line fitting are available in the Ap- 
pendix. The methods used to develop and apply the template are also described. We 
illustrate the application of the derived Fe n and Fe in templates by fitting and subtract- 
ing iron emission from the spectra of four high redshift quasars and of the nearby quasar, 
3C273, confirming their general applicability to active galaxies despite the somewhat 
unusual properties of IZwl. We briefly discuss the small discrepancies between the 
observed iron emission of these quasars and the ultraviolet template, and compare the 
template with previously published ones. We discuss the advantages and limitations of 
the UV Fe n and Fe in templates and of the template fitting method. We conclude that 
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the templates work sufficiently well to be a valuable and important tool for eliminating 
and studying the iron emission in active galaxies, at least until accurate theoretical iron 
emission models are developed. 

The Si iv+O iv] A 1400 feature in I Zw 1 is clearly strong relative to C iv A 1549, and 
Civ and Cm] A 1909 are both relatively weak. This may partially be due to the higher 
densities and lower ionization parameter prevailing in Narrow Line Seyfert 1 galaxies, 
and to the big blue bump shifting towards lower energies in more luminous Seyferts, 
such as I Zw 1 . In I Zw 1 the narrow line width reveals that C m] is heavily blended with 
Si III] A 1892, Aim AA 1854,1863, and Fein transitions. This suggests that the Cm] line 
strength and width may be overestimated in many quasar line studies where the lines 
are broader and deblending is not possible. This affects density estimates of the broad 
line region. Photoionization modeling, including all these line features, and subsequent 
fitting to the spectra are required to estimate the true C ill] strength. We also argue, 
based on earlier work, that (strong) iron emission may be connected with high densities 
and associated with outflows. 

Subject headings: galaxies: active — galaxies: Seyferts: individual (I Zw 1) — methods: 
data analysis — quasars: emission lines 



1. Introduction 

Quasar 5 ultraviolet (hereafter UV) spectra characteristically contain broad emission lines orig- 
inating in the central (< light year) region, the broad line region (BLR). The strongest lines are (the 
resonance lines) Lya A 1216, Si iv+O iv]A 1400, CivAA1548, 1551, Cm] A 1909, and MgnAA2796, 
2803. For a number of years after the discovery of quasars these lines and a few iron transitions 
(Greenstein & Schmidt 1964; Wampler & Oke 1967) were the only lines detected in their UV 
spectra. With the availability of high signal-to- noise ratio (S/N), high resolution spectra, and in 
particular data taken with Hopkins Ultraviolet Telescope (e. g., Kriss et al. 1992; Zheng et al. 1995; 
Zheng, Kriss, & Davidsen 1996), Hubble Space Telescope (HST) (e. g., Laor et al. 1994, 1995; Zheng 
et al. 1997; Kriss et al. 2000; Kraemer &: Crenshaw 2000), and the Keck telescope (e. g., Tran, Co- 
hen & Goodrich 1995; Brotherton et al. 1997; Barlow & Sargent 1997; Larkin et al. 2000; Carson 
et al. 2000) it has become obvious that quasar spectra contain a plethora of weak lines as well. 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data archive at the 
Space Telescope Science Institute. STScI is operated by the Association of Universities for Research in Astronomy, 
Inc. under the NASA contract NAS 5-26555. 

5 "Quasars" here refers to both the radio-loud and radio-quiet subgroups 
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1.1. The Need for Fitting and Removal of Iron Emission 

In order to study e. g., the relationship between the continuum and the line emission, the 
dynamical structure, the ionization balance and structure, or the chemical abundances in active 
galactic nuclei (herafter AGNs) and quasars it is important to make reliable measurements of the 
continuum emission, the strong broad lines including the wings of their profiles, and the weak 
emission lines. In addition to hydrogen, helium and the above-mentioned elements, AGN and 
quasars also contain atomic and ionic iron, the stable end product of nucleosynthesis. Due to the 
large number of electron levels in iron atoms, thousands of emission line transitions are distributed 
throughout the UV and optical spectral regions. The weak lines, iron as well as those not being iron 
transitions (hereafter 'non-iron' lines), blend together, partly due to the associated transitions being 
very close or overlapping in energy, and partly due to the high (presumed) dynamic velocities of the 
broad-line emitting clouds broadening the lines. Line widths [F WHM, the full width at half peak 
flux] of up to 10 000 km s _1 have been measured (Wilkes 2000). This heavy blending of weaker lines, 
dominated by iron transitions, forms a pseudo-continuum above the intrinsically emitted continuum 
(see e. g., Figs. 6 and 7) even in high S/N quasar spectra (Wills et al. 1985; Boroson & Green 1992, 
hereafter BG92; Wills et al. 1995). This pseudo-continuum severely complicates the study of both 
weak non-iron features, particularly important for abundance and ionization studies, and the wings 
of strong resonance lines. The uncertainty introduced into broad emission-line measurements by 
iron emission "contamination" can be relatively large: excluding very weak lines where it can 
reach ~100%, a rough estimate is 5—50% based on Wilkes (1984) and our own line measurements 
before and after iron emission removal using the template presented here (§ 4). The exact level of 
uncertainty depends on the line transition, the line parameter measured and the overall strength of 
the emission line spectrum. Moreover, the relative strengths of the iron features vary greatly from 
object to object. These uncertainties are dominated by the difficulty in determining an accurate 
continuum level, one of the main problems in quasar line studies, although such errors are seldom 
quoted. Combined, these uncertainties provide a strong argument that all quasar spectral studies 
should include iron emission fitting and removal. 

These limitations imposed by the contaminating iron emission have been known for some 
time (e. g., Wills Sz Browne 1986). A possible solution to the continuum level uncertainty is to 
use only wavelength ranges which contain pure continuum emission (so-called continuum windows; 
Francis et al. 1991) when fitting the continuum. However, few such regions exist and they are 
generally small, especially in the UV spectral region, ~ 1000 — 2000A. Moreover, in objects with 
very broad lines even these line-free regions may be contaminated by broad line-wing emission 
(Boroson, Persson & Oke 1985; Wills & Browne 1986). 
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1.2. Earlier Work on Iron Emission 

Allowance for the optical and UV iron emission has not extensively been made, inspite of the 
increasing recognition of its importance for broad line studies. It was rendered a very hard task in 
the 1980's by the lack of sufficient atomic data to allow identification of the iron transitions and 
high spectral resolution to resolve them. Earlier work includes identifying the UV iron emission 
lines in AGN [e. g., Wills, Netzer k Wills 1980; Penston 1980; Penston et al. 1983 (in stars); 
Hartig k Baldwin 1986, Johansson k Jordan 1984] and modeling this emission (e. g., Netzer 1980; 
Kwan k Krolik 1981; Netzer k Wills 1983; Wills et al. 1985; Collin-Souffrin et al. 1986; Penston 
1987; Collin-Souffrin, Hameury k Joly 1988; Krolik k Kallman 1988; Ferland k Persson 1989; 
Dumont k Collin-Souffrin 1990; Netzer 1990, and references therein). Despite extensive efforts the 
current theoretical models of the optical and UV iron emission cannot fully reproduce and explain 
the observations. Fortunately, projects are under way to compute detailed radiative transition 
probabilities of iron atoms and ions (the IRON Project: Hummer et al. 1993; Nahar, Bautista,& 
Pradhan 1997; Nahar et al. 2000, and refererences therein) and compute photoionization models 
of AGN which also take the iron emission into account with or without exact radiative transfer 
(Sigut k Pradhan 1998; "The Kentucky group": Verner et al. 1999, and references therein) based 
on improved atomic data (e. g., the Opacity Project: Seaton et al. 1994). Other studies presenting 
iron line lists and radiative transition probabilities, which are also used here, include Fuhr, Martin 
k Wiese (1988), Ekberg (1993), Giridhar k Arellano Ferro (1995), Nahar (1995), Nahar k Pradhan 
(1996), Quinet (1996), Quinet, Le Dourneuf, k Zeippen (1996), and Kurucz k Bell (1995). These 
iron line lists are available on the world wide web. 

1.3. Iron Emission Correction 

Phillips (1977) was one of the first to compare the spectrum of IZwl to other Seyfert galaxy 
spectra by broadening the former with a Gaussian profile, selected to match the line widths in the 
latter. Phillips did not subtract the iron emission but confirmed his suspicion that the Fe n and H I 
lines have essentially the same widths and profiles. Wills et al. (1985) simulated the Fen emission 
and absorption in eight low to intermediate redshift quasars using photoionization modeling of 
~3000 iron lines. In spite of the limitations imposed by the atomic data and computing speed at 
the time, they were able to reproduce the UV-optical spectra reasonably well. Some deviations 
from the observed spectra are, however, seen (see § 3.4). The photoionization models presented by 
Wills et al. are specific to their individual object spectra and were not generated with a general 
subtraction of iron emission in AGNs in mind, where a typical iron spectrum is more appropriate. 
Boroson, Persson k Oke (1985) made one of the first attempts to take the effects of the iron emission 
on broad- line studies into direct account. They estimated the uncertainty associated with line 
(equivalent width, EW) measurements as a function of line widths by broadening an AGN spectrum 
by a Hanning profile (triangular) with a range of widths. They found the measured EWs to weaken 
with increasing FWHM because the continuum level becomes systematically overestimated, as 
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the iron lines broaden and blend to form a pseudo-continuum. As a consequence, narrow lines 
incorrectly appear stronger. This may explain why quasars, with their broad emission lines were 
not recognized early on to be strong iron emitters 6 (Phillips 1977; Davidson & Netzer 1979). One 
way of correcting for the iron emission is to model it using the profile information of the H (3 line, 
as the isolated Fe n and H (3 profiles are observed to have both similar widths and profile shapes 
(Phillips 1977; Laor et al. 1997b; hereafter L97). Thus a synthetic spectrum of the iron emission 
can be constructed by shifting and scaling H (3 profile templates according to a list of iron line 
positions and relative strengths. However, this method depends highly on the atomic data lists to 
be representative, in terms of accuracy and completeness, of the iron emission in AGN and quasars. 
As mentioned, theoretical studies are not yet able to fully reproduce the observed iron emission, 
but work is in progress (D. Verner 1998, private communication; Verner et al. 1999; see also § 3.4). 

A significant improvement in the accuracy of broad emission line measurements can be ob- 
tained by using an observed AGN 7 or quasar iron emission template to fit and subtract the iron 
emission in quasar spectra before performing line and continuum measurements. The benefits of 
using such a template, containing all the iron transitions typically present in quasars, are manifold. 
Combined with the much improved atomic data, becoming increasingly available (§ 1.2), empirical 
templates provide a powerful tool to study the observed iron spectrum in terms of theoretical mod- 
els. This template method was adopted by BG92, who successfully used an optical (4250— 7000A) 
iron template based on the nearby, narrow-line Seyfert 1 (NLS1; Osterbrock & Pogge 1985), IZwl, 
on the Bright Quasar Survey sample. Their method is described in § 4.2. Corbin & Boroson (1996) 
use the same method with a 2300-3000A UV iron template of IZwl. See § 3.4 for a comparison 
with this template. 

The advent of HST UV quasar spectral data allows us to extend the template method into the 
UV regime. We here present, to our knowledge, the first high S/N, high resolution, quasar empirical 
UV iron template spectrum ranging from restframe AA1250 to 3090A which is applicable to quasar 
data. The template is based on HST (archival) data of IZwl (L97). The method with which the 
template was generated is an extension of that of BG92 and Corbin & Boroson (1996) including a 
more detailed fitting of the lines which are not iron (see § 3 for details and § 3.4 for a discussion). 
When optical data between 3000A and 4250A become available (L97 present 3000— 3800A data), 
suitable for production of an intermediate wavelength iron template, the available templates will 
provide a powerful tool for consistent modeling of the UV and optical iron emission throughout the 
1250A to 7000A region of quasar spectra. 



The fact that radio-loud quasars were often studied more than radio-quiets, due to their easily-detectable, strong 
radio emission, may also have contributed, because radio-quiet quasars show relatively stronger optical iron emission 
(Peterson, Foltz & Byard 1981; Bergeron & Kunth 1984; Corbin 1997). 

7 The formal distinction between the use of the term 'AGN' and 'quasar' is only that of the object's luminosity, 
with quasars occupying the most powerful end of the luminosity range of active galactic nuclei (Mv < —23, Veron- 
Cetty & Veron 1993). As the methods and iron templates described here are applicable to AGNs as well, this fact is 
silently assumed in the following, when we refer to quasars only, in addressing the use of the template. 
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The procedure of fitting the iron emission using an empirical template assumes that (1) the 
iron spectrum of I Zw 1 is representative of quasar iron emission spectra and that (2) all quasars 
have similar iron spectra to within a scaling factor and/or a line profile broadening. Not enough 
is currently known about the iron emission in AGN to firmly assess the validity of these assump- 
tions. Though until theoretical models can explain the observed iron emission more confidently this 
template fitting procedure is the best available approach. We test the method by fitting several 
representative quasar spectra (§ 4.3) showing that the iron emission can be well fitted if allowed the 
freedom to vary some of the multiplet ratios in the empirical template, indicating that the basic 
assumptions are valid. 

The iron emission templates have importance not only for our ability to fit and subtract the iron 
emission in quasar spectra, but also as tools with which we can study the iron emission strengths 
themselves. Iron is a key coolant emitting ~ 25% of the total energy output from the BLR (Wills 
et al. 1985; Boiler, Boiler, Brandt & Fink 1996), emphasizing the importance of including the iron 
emission in studies of the BLR. 

1.4. The Narrow-Line Seyfert 1, I Zw 1 

The spectra of NLS1 galaxies are particularly useful for generating empirical iron templates 
because their strong and rich iron emission allows detection of weak iron features and identification 
of as many iron transitions (Figs. 2 and 4, Table 4, § 3.3.1) as are typically present in AGN 
spectra. The relatively narrow width (FWHM<2000kms^ 1 ) of the broad emission lines permits 
the individual non-iron lines to be resolved, isolated, and removed from the spectrum. It also allows 
us to match the line width in most iron-contaminated AGN spectra so as to fit and subtract the 
iron emission by broadening the iron template. 

IZwl (PG 0050+124; z = 0.061, see § 2) is classified as a NLS1 and has the "narrow" broad 
emission-lines, strong Fen emission (e. g., Phillips 1976; Osterbrock & Pogge 1985; Lipari, Ter- 
levich, & Macchetto 1993; Pogge 2000; Rodrfguez-Ardila, Pastoriza & Donzelli 2000; Rudy et al. 
2000; but see also Gaskell 2000), steep soft X-ray spectra (Boiler et al. 1996), and strong far-IR 
emission (Halpern &: Oke 1987), typical of this class. Their radio properties, however, are similar 
to those of other Seyfert galaxies (Ulvestad, Antonucci & Goodrich 1995). They are also variable 
at optical and UV wavelengths (e. g., Zwicky 1971; Giannuzzo & Stirpe 1996; Rodriguez-Pascual, 
Mas-Hesse, & Santos-Lleo 1997; Leighly 1999; Miller et al. 2000) and the fastest X-ray variable 
AGN known (see discussion by Boiler et al. 1996). I Zw 1 is a well studied object thanks in part to 
its relative brightness (Schmidt & Green 1983), its exceptionally narrow lines, its iron-rich spectrum 
(Sargent 1968; Phillips 1976, 1977; Oke & Lauer 1979; Boroson et al. 1985; Halpern & Oke 1987; 
BG92; L97), and its richness in low ionization lines (Persson & McGregor 1985; van Groningen 
1993; L97; Table 3). IZwl is also an infrared (IR) luminous source (Rieke Sz Low 1972; Rieke 
1978; Halpern & Oke 1987) and has been observed in CO and H2 molecular lines (Barvainis, Alloin 
& Antonucci 1989; Eckart et al. 1994). Some modeling of the observed optical iron emission was 
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performed by Phillips (1978). IZwl was observed by IUE in 1978-1982 (Courvoisier & Paltani 
1992; Lanzetta, Turnshek & Sandoval 1993; Paltani & Courvoisier 1994), and the source is included 
in a number of IUE studies addressing larger samples of AGN (e. g., Wu, Boggess &; Gull 1983; 
Pian & Treves 1993; Wang, Zhou & Gao 1996; Rodn'guez-Pascual et al. 1997). Condon, Hutchings 
& Gower (1985) detected 21 cm radio emission from I Zw 1 consistent with emission from a late-type 
host galaxy. 

IZwl is a particularly good choice for an empirical Fen template as it is so well studied, 
especially in terms of its optical iron emission and was used by BG92 and Corbin & Boroson (1996) 
for their templates. 

Additional candidates of narrow-line, iron-rich AGN suitable for use as iron templates may 
be found among other NLSls. Possibilities include Mrk957 (5C 3.100), Ark564, IE 1226.9+1336, 
E1228+123, Mrk507 (1748+687), Mrk42 (1151+465), and the less extreme Mrk291 (1552+193), 
Mrk493 (1557+352) and 1244+026. 

The structure of this paper is as follows: § 2 addresses the data processing, § 3 contains a description 
of the generation of the template and identification of the various spectral lines, in § 3.4 comparisons 
are made with other available empirical UV templates and synthetic Fe n models, § 4 describes the 
application of the template along with initial results, § 5 is dedicated to a discussion of I Zw 1 and 
some of its spectral features, and § 6 summarizes the main conclusions. Comments on individual 
line features are deferred to Appendix A. 

2. Data Processing 

The data presented here consist of HST Faint Object Spectrograph (FOS) archival spectra of 
I Zw 1. The journal of observation and instrumental setup is summarized in Table 1 for convenience 
(see also L97). 

The G130H spectrum was calibrated according to the standard CALFOS procedure with up- 
dated (1996 March) calibration files, as the flux calibration status of the archival data is uncertain 8 . 
The pipeline calibrations were used for the G190H and G270H spectra. Multiple observations ob- 
tained with the same grating were combined by weighting with the exposure time to form a single 
spectrum for each wavelength region. A color excess, E(B— V) = A#/48 x 10 20 cm -2 = 0.105 mag, 
was determined based on the Galactic hydrogen column density, Nh = 5.05 (±0.1) x 10 20 cm~ 2 , 
observed 9 by Stark et al. (1992). The spectra were dereddened using the average extinction curve 



8 The data were skipped by the pipeline calibration, evident from the data intensity level and the missing flux 
calibration flags in the data headers. 

9 Note: Errors are not quoted by Stark et al. (1992), but the applied correction technique for stray radia- 
tion was developed by Lockman, Jahoda & McCammon (1986), who quote an error of ~1 xl0 19 cm~ 2 . Using 
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presented by Cardelli, Clayton k Mathis (1989), using Ay = 3.1 * E(B-V) and the IRAF 10 (V2.ll) 
task 'deredden'. Due to the average nature of this curve, some residual dust extinction features 
may be present in the spectra, especially in the 2200A region (Fig. 1), where characteristic dust 
extinction features are generally expected (e. g., Osterbrock 1989). Iron emission in this region com- 
plicates this identification. We corrected for an offset in the wavelength solution, due to the non-zero 
uncertainty in the wavelength calibration (typical uncertainties are ~0.25A in G130H, ~0.37A in 
G190H, ~0.52A in G270H; Leitherer 1995), by comparing the observed wavelengths of the galac- 
tic interstellar medium (ISM) absorption lines with their laboratory wavelengths. The absorption 
lines used are: Sin A1190, Sin A1193, Sim A1206, Sin A1260, CnA1335, Sin A1527, AlnA1671, 
FenA2344, FeiiA2374, FeiiA2382, FeiiA2586, FeiiA2600, Mgn AA2796, 2803 and MgiA2853 (Ta- 
ble 2). The applied wavelength offsets are listed in Table 1. These offsets, consistent with those 
applied by L97, result in the absorption line positions matching the laboratory wavelengths (e. g., 
Savage et al. 1993; Morton 1991) to within ±0.3A (rms). 

The spectra from the individual gratings were co-added to produce a full 1141 A - 3278A 
(1075A - 3090A restframe) spectrum. The spectra from gratings G190H and G270H were obtained 
the same day and show no difference in the continuum level, and so were co-added without scaling 
either spectrum, using the average flux in the overlapping region. The continuum-level in the G130H 
spectrum, taken 6 months earlier, is clearly shifted relative to that of G190H. The G130H spectrum 
normalization is described in § 2.1. For consistency, the combined, final spectrum was rebinned 11 to 
match the dispersion of the G270H grating data, AA = 0.51lA/pix, the lowest dispersion available 
in the three gratings (Table 1). The spectrum ranges from 1075.2 to 3089. 8A, has a dispersion of 
0.482 A/pix and a resolution of 1.86A in the restframe (Leitherer 1995). The IZwl spectrum is 
shown in L97 and Figure 1 (dotted line). A redshift of z = 0.061, defined by the Mgn profile peak, 
is used throughout this paper, and is consistent with measurements by Phillips (1976), Schmidt k 
Green (1983), Wu et al. (1983), Persson k McGregor (1985), van Groningen (1993), Condon et al. 
(1985) and Smith et al. (1997). Condon et al. determine a redshift of 0.061136 to an accuracy of 
8xl0 -6 using Hi radio measurements. 

The spectral energy distributions of quasars and AGNs (Elvis et al. 1994) span 11 orders of 
magnitude in frequency from X-rays to millimeter wavelengths (and to radio wavelengths for some 
objects; Weedman 1986; Peterson 1997). Observations (e. g., Oke, Shields k Korycansky 1984; 



the same correction technique Elvis, Lockman & Wilkes (1989) measure an HI column density towards IZwl of 
Nh = 5.07 (±0.1) x 10 20 cm~ 2 consistent with the Stark et al. measurement. The slight offset between the two Nh 
measurements will not significantly affect the reddening correction. 

10 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science 
Foundation. 

"The rebinning is necessary in order to combine and simultaneously process the G130H, G190H, and G270H 
spectra for generating the template. The resultant resolution is sufficient for the purpose of generating and applying 
the template. 
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Wills et al. 1985) suggest that the continuum at optical and UV wavelengths can be approximated 
by a single power-law, F„ ~ u~ av . We thus chose to fit a power-law continuum to the IZwl 
restframe spectrum before any emission line features were fitted. After the initial completion of 
the iron template, which included careful fittings (§ 3.3.2) of the emission lines which are not 
iron, it was realized that due to an STSDAS software bug 12 the archival G130H spectrum was 
unknowingly calibrated with the wrong (pre-costar) flux calibration files. These data were then 
recalibrated 13 using the most recent (1996 March) calibration files. We then fitted a power-law 
continuum from 1075 to ~1720A, independent of the power-law continuum at longer wavelengths 14 . 
The resulting continuum is a broken power-law with the break at 1716A: a blue continuum slope, 
a v = 1.9 and normalization, F^(1500 A) = 3.45xl0 -14 ergs cm~ 2 s _1 A -1 , and a red continuum 
slope, a„ = 1.0 and F A (1500A) = 3.89xl(T 14 ergs cm 2 s 1 A 1 . The choice of this continuum 
(Fig. 1) does not affect the application of the template, as it will be scaled and broadened to match 
the target spectrum in the application process (§ 4). The continuum windows used in the fitting 
(AA1312-1327, 1347-1353, 1641-1647, 1675-1690, 3007-3027) are specific to the I Zw 1 spectrum 
and are different from those suggested by Francis et al. (1991), which were based on an average 
AGN spectrum. The spectrum of I Zw 1 is rich in low-ionization lines and weak emission features 
(Fig. 1; L97) which contaminate the "average" continuum windows. 



2.1. The G130H Spectrum Scaling and IZwl Variability 

The variation in continuum level between the G130H and the G190H+G270H spectra taken 
six months apart is consistent with the known variability of IZwl (e. g., Zwicky 1971; Giannuzzo 
& Stirpe 1996; Rodrfguez-Pascual et al. 1997; Leighly 1999). The normalization of the recalibrated 
G130H spectrum to the G190H+G270H spectrum was determined as the ratio of the median flux in 
the overlapping region of the spectra (1.3±0.56). As the bluest ~40 A of the G190H spectrum shows 
a relative error of 42% as opposed to 4% in G130H, the normalized G130H spectrum substituted 
that of G190H in this overlapping region. However, the G190H flux uncertainty entirely dominates 
the formal error on the normalization, and so useful line flux estimates cannot be deduced below 
1500 A. In spite of the large formal error, a smooth spectrum and power-law continuum resulted 
in the A1075-1716A region, implying that the actual uncertainty is lower. As the template is 
best applied by subdividing the A1075-3090A range before scaling to match the iron emission 
strength in individual quasar spectra (§§ 4.2 - 4.3), its use is not adversely affected by the G130H 
normalization. In fact, the poor match of the Fe II emission in the G130H grating spectrum for some 

12 the relevant task, "ADDNEWKEYS" , has now been updated by the STScI HST help desk. 

13 This recalibration results in a flux level ~9 % higher than that of L97, because the FOS calibration files were 
updated (1996 March) later than their data processing in 1994. 

14 The emission and absorption models below 1700A were then regenerated and new iron templates (of Fell and 
Fein emission) were constructed. 
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objects (see § 4.3) directly shows the need to subdivide the template in order to obtain optimum 
fits to AGN spectra. This is because all the individual iron multiplets did not brighten with the 
same factor in IZwl as did the continuum in the time between the G130H and the G190H spectra 
were observed. If iron multiplet strengths vary in the same object, they are also very likely to vary 
among objects. This is also evident from the fact that not all AGNs have Fen UV 191, Fein UV34, 
and/or Fein UV47 multiplets as strong as does IZwl (§ 4.3 and Appendix A). Lanzetta et al. 
(1993) present a spectrum of IZw 1 representing the average of the UV spectra observed with the 
International Ultraviolet Explorer (IUE) from 1978 through 1982. Comparison with the HST data 
to confirm the scaling of the G130H spectrum is, however, not possible due to the widely differing 
host galaxy contribution 15 and to line and continuum variations (cf. Wu et al. 1983; Pian & Treves 
1993; Wang et al. 1996; Rodriguez-Pascual et al. 1997; L97). 

The only other UV data available are, to our knowledge, spectropolarimetric data in the HST 
data archives. A comparison with such data requires knowledge of the scattering medium and an 
understanding of the scattered spectrum, which are beyond the scope of this work. Hence, it is not 
possible at present to constrain the absolute scaling of G130H spectrum further, though as noted 
this uncertainty does not adversely affect the applicability of the iron template. 



3. The Iron Emission Template 

3.1. Development of the Template 

Once the data were calibrated, rebinned to a common dispersion, and co-added the template 
was generated using the following procedure: 

• A power-law continuum was fitted to pure continuum wavelength regions (§ 2) in the spectrum 
and this continuum fit was subtracted (result is overplotted in Fig. 1; dotted line) 

• All absorption and (non-iron and iron) emission features were identified, including the strong, 
weak, blended and unblended ones (Fig. 1, Tables 2 and 3; § 3.2, § 3.3.1) 

• The non-iron emission features (Fig. 5) and Galactic/ISM absorption features (§ 3.2, § 3.3.2)were 
fitted. This required a simultaneous fitting of some Fein features (Table 4). Separate ab- 
sorption and emission models were created (Fig. 1) 

• The remaining Fein features were fitted and deblended, as needed (Table 4). The Fein line 
emission was isolated and an Fein model was created (§ 3.3.2) 

• The absorption and non-iron emission models were subtracted from the original spectrum to 
create an iron template (containing both Fe n and Fe in emission; Fig. 2) 



15 The IUE spectrograph has a large aperture (10" x 20", Boggess et al. 1978) thereby sampling the AGN host 
galaxy. 
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• The Fe in emission model was subtracted to create a pure Fe n template (Fig. 3) 

• The pixel values were set to zero in the two iron templates in regions containing pure noise 
residuals left over from subtracting emission and absorption feature fits (Fig. 2 and 3). This 
prevents introduction of artifacts and noise to the target spectra when the template is applied. 

In the following sections we discuss the fitting and identification steps of this procedure in 
more detail. We also compare the template, in § 3.4, with other UV iron templates and theoretical 
models, currently available. 

3.2. Absorption Lines 

Features were identified as absorption lines (Table 2) when the minimum flux in the feature 
deviated by more than 3 cr (a = rms of the fluxes around the average local spectrum level) from the 
continuum or if the position of the feature coincided with that expected as part of a doublet line 
where the strongest line component was already identified. The identification of the specific ion and 
transition responsible for each absorption line was made using the list of typical ISM absorption 
lines by Savage et al. (1993) and the line lists by Morton (1991). All the absorption lines detected 
and identified in the spectrum are consistent with Galactic absorption with the exception of the 
features at 1306A and 1310A which are due to N v AA 1238, 1243 (blueshifted) absorption associated 
with IZwl (L97). 

Absorption features were fitted by multiple Gaussian components assuming a constant width 
for all components of a multiplet. The fit was then subtracted from the spectrum. 

3.3. Emission Lines 

3.3.1. Line Identifications 

Our identifications of each emission feature are based on reference line lists for non-iron (Wilkes 
2000, Morton 1991, Verner, Barthel, & Tytler 1994), and iron transitions [Moore 1950 (Fell and 
Fein multiplets); Penston et al. 1983 (Fen); Fuhr et al. 1988 (Fell and Fein); Giridhar & Arellano 
Ferro 1995 (Fen); Nahar 1995 (Fen); Nahar & Pradhan 1996 (Fein); Kurucz & Bell 1995 (Fen 
and Fein); D. Verner (1996, private communication, Fen)]. L97 suggested identification of a 
large number of features in the I Zw 1 spectrum, but that paper was not our main reference as it 
appeared after our work was commenced and after most of the line identifications were completed. 
The data from L97 are essentially the same data presented here. Slight differences exist as this 
G130H spectrum is recalibrated with more recent files and is renormalized to the level of the 
G190H spectrum (see § 2). Certain differences are present between our work and that by L97: 
(1) we identify a few additional non-iron features, (2) we suggest identifications of individual iron 
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transitions in addition to the multiplet identifications in L97, (3) we find slightly different velocity 
shifts of the various line groups, and (4) many of the L97 line measurements deviate, though by no 
more than 30% for the stronger and/or isolated lines. We briefly discuss (3) and (4) in § 5.1.2 and 
§ 5.1.3. 

In Appendix A we comment on the individual line features and their fits. A thorough discussion 
of line intensities is given by L97 and not repeated here apart from brief discussions of the weak 
Cm] (§ 5.2) and the relatively strong A1400 emission (§ 5.3). 

Non-Iron Emission Lines The UV HST spectrum of IZwl has a sufficiently high S/N and 
spectral resolution to permit detection and identification of many weak features in addition to the 
strong, broad emission-lines commonly observed in quasars. Table 3 lists the detected (non-iron) 
line features and their identification along with basic line parameters; see § 3.3.2 on line fitting for 
further details. L97 also identify most of these lines. We contribute with a few extra identifications 
and measurements: Sim* A1297, and Sin AA 1527,1533, based on Verner et al. (1994), Morton 
(1991), and Wilkes (2000). We are not able to confirm the [Cm] A1907 feature identified by L97 
due to our slightly degraded resolution (§ 2). 

Iron Emission Features Figure 4 shows suggested identifications of individual, mostly un- 
blended Fen and Fein emission UV multiplets 16 . The identifications are based on a visual in- 
spection of clearly visible peaks in the spectrum whose positions and relative strengths were cross- 
correlated with the multiplet table of Moore (1950). This goes a step further than the work by L97, 
who mark in their figure 2 only the expected wavelengths of some of the strongest iron multiplets. 
The length of the marker of a given transition in a given multiplet is proportional to its oscillator 
strength (the 'intensity' listed by Moore is a rough measure of the relative oscillator strengths in 
the multiplet; note, that the physical conditions folded with the oscillator strengths will give the 
observed line strengths and ratios). The scaling factor is the same for all multiplets (and is arbi- 
trarily chosen). Note that the transitions in the spectrum sometimes appear blueshifted by 1-2 A 
relative to the laboratory wavelengths. See e. g., FemUV34 at ~1914A. Each label contains first 
the ionization level, then the UV multiplet number, separated by a hyphen. That is, '2-104' denotes 
the Fen UV104 multiplet, while '3-158' denotes the Fein UV158 multiplet, etc. Fen and Fem 
multiplets in the spectrum for which the relative transition strengths in the multiplet do not appear 
to follow the ~optically thin multiplet strengths listed in Moore (1950) have an 'm' attached to the 
labeled multiplet number. A '?' indicate that this multiplets presence is uncertain. Labels of mul- 
tiplets which appear slightly blended are shown in parentheses. Square brackets denote multiplets 
whose presence is suggested by the fitting process (see Appendix A); Fen UV10, which coincide 
with the geocoronal Lya emission, is expected only. Heavily blended multiplets are not labeled, 
especially in the small blue bump region from ~2650 A to 3090 A. The lack of labeled multiplets 



A full scale version of Fig. 4 can be found at http://www.astronomy.ohio-state.edu/~vester/IronEmission. 
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between Ly a and C iv does not denote a lack of Fe n or Fe in emission transitions. Most transitions 
in this range are too weak (relatively) to have multiplet numbers assigned. Blueward of Lya an 
accurate identification of iron emission multiplets requires advanced modeling of the spectrum due 
to the heavy blending with Ly a forest and other absorption lines. This is beyond the scope of this 
paper. L97 mark possible iron multiplets in this region. No obvious Fe II or Fe III transitions arc 
identified immediately redward of Civ (~1550 - 1700 A). Marziani et al. (1996) discuss some Fell 
multiplets in this wavelength region in earlier HST data of I Zw 1 (their figure 2). Some singlet Fe n 
and Fe in emission features are identified and commented on in Appendix A. 

3.3.2. Fitting Procedures 

The emission lines, listed in Tables 3 and 4, were fitted using single or multiple Gaussian 
components, as needed, and then subtracted from the spectrum. The Fein and Fen features were 
separated by use of the fitted Fe m lines in Table 4 to generate two templates, one for each ion 
(§ 3.1). The IRAF task 'splot' was used for the Gaussian fitting due to its interactive nature and 
the ease with which the fitted components can immediately be extracted and compared to the 
data, 'splot' is often thought of as a tool for first order estimates but its characteristics allow us to 
constrain our fits faster and more conveniently than other immediately available fitting programs, 
due to their largely non-interactive nature. 

Fitting the multiplets and line complexes was an iterative process. First the narrow components 
were fitted to the visible part of the lines and subtracted, then the broad components were fitted to 
the residuals. The broad component fits were then subtracted from the original data permitting an 
improved fit to be obtained for the narrow lines, which in turn were subtracted to improve the broad 
component fit, thereby iterating to obtain an optimal solution. Similar iterations were performed 
for individual (narrow) components in doublets and in regions where the narrow components are 
somewhat blended (e. g., Civ, Mgn, and the Cm] complex). When a line complex contains two 
broad components (the Aim and Mgn-doublet complexes, Figs. 5d and 5g), they are often too 
heavily blended to be well separated and so were fitted by a single broad component. Two broad 
components were fitted simultaneously to Sim] and Cm], but as they are heavily blended the fit 
is not unique (see below; Fig. 5d; Appendix A). 

The average line width of the singlet lines is 900 ±150 km s _1 (Tables 3 and 4) consistent with 
previous studies of the line emission in IZwl (Phillips 1976; BG92; van Groningen 1993; L97). 
However, lines with widths as small as ~ 300 - 400 kms -1 (the spectral resolution limit) are 
detected, as discussed in § 5.1.1. These are often Fen and Fein lines (see Table 4). 

Each Gaussian component is defined by three parameters: position, width, and strength. In 
the heavily blended line complexes the absolute strengths of the individual Gaussian components 
are not well constrained. The uncertainty in the intensity of individual lines is estimated at 10—50% 
depending on the width and strength of the line and on how well the data constrain the fit. The 
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more blended the component, the more uncertain is the fit. Due to the large number of parameters 
involved in fitting the large line complexes, individual component solutions are not unique. In 
such cases, those solutions yielding parameters consistent with other lines were preferred. We 
note that, when the strength of any line fit was not well constrained by the observed 
spectrum, the affected iron residual emission was purposefully underestimated in order 
to avoid overcorrection of the quasar spectrum to which the template is applied. The 
representation of the entire line complex (i. e., the sum of all the individual components) was, 
however, always well constrained by the data and the estimated errors are of order a few percent. 

Examples of the emission- line fits can be found in Figure 5. The individual measurements of 
fitted non-iron and iron lines are listed in Tables 3 and 4, respectively. Measurements are based 
on the individual Gaussian components making up the line profile; sums for the line complexes are 
also listed. The equivalent widths are always measured relative to the adopted (global) continuum 
level (§ 2). Very weak iron features were not fitted. 

We emphasize that the multiple Gaussian component fitting is simply a tool, and we make 
no assumptions as to the mechanism responsible for the overall profile shapes. Thus, as no phys- 
ical meaning is associated with the individual Gaussian component fits, the absolute and relative 
fluxes of the fitted components cannot be interpreted in terms of physically distinct emitting re- 
gions. Extreme care should always be exercised when interpreting the measurements, particularly 
of weak features. Features for which certain fitting parameters and/or the line identifications are 
particularly uncertain are mentioned in Appendix A. 

3.4. Comparison with Earlier UV Iron Templates and Models 

Despite the wealth of work on modeling the iron emission in quasars there is no UV template 
electronically available which can be applied to observed quasar spectra, including our own (M. 
Vestergaard et al. 2000, in preparation), for removal of the UV Fen emission. One other empirical 
UV template is available in the literature: the Corbin & Boroson (1996) template covering the 
iron bump around Mgll. Below we briefly comment on how this compares to our template in the 
overlapping region. A couple of synthetic Fell spectra have been generated based on the available 
knowledge of the iron emission mechanism, the quality of the atomic data, and the computing facil- 
ities at the time (Wills et al. 1985; Verner et al. 1999). We briefly discuss the apparent differences 
with the current iron template. 

3.4-1- Empirical UV templates 

Corbin & Boroson (1996) present a AA2300-3000 A iron template also derived from the HST 
spectrum of IZwl (L97). Although a detailed comparison is not possible (the digital spectrum 
is not available and Figure 2 in Corbin & Boroson (1996) shows an already broadened version of 
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the template), it is clear that differences are present due to their use of interpolation rather than 
deblending to remove non-iron features. This is particularly clear around Mg n where interpolation 
across the Mgn line likely overestimates the iron emission strength by ~300%; the Fen emission 
level is at most 15-20% of the peak height of the Fen on either side of the Mgn gap (cf. figure 12 
in Verner et al. 1999). In contrast, our approach (§ 3.3.2) is to purposefully underestimate the iron 
emission when the strength of the iron and/or the non-iron emission is in doubt. We have therefore 
assumed that the broad component of Mgn, the existence of which is evident in Figures. 1 and 5g, 
dominates this wavelength region, consistent with the findings of L97, and so our template shows 
a gap in the iron emission around Mgll. 

Our sample fitting and subtraction (§ 4.3) of the iron emission in spectra of high-redshift 
quasars using the template presented here yields residual profiles of the Mg n line resembling those 
of the other prominent UV lines, supporting our method. Our Gaussian fitting of Mgn indicates 
that one will underestimate its strength by a factor of ~2 if the broad Mgn component is ignored, 
as in the interpolation process. 

3.4-2. Synthetic Fell Spectra 

Wills et al. (1985) present theoretical models of the 1800-5000 A Fen emission, the Balmer 
continuum and lines, and compare with observed spectra of a small sample of low-redshift quasars 
(z ~0.12-0.6). Their success in simulating the observed spectra is generally good, though a few 
discrepancies in line strengths among multiplets remain. Over the past few years the Opacity and 
IRON projects have applied the power of modern computers to determining the atomic parameters 
of the thousands of possible iron transitions (e. g., Hummer et al. 1993, Seaton et al. 1994, Nahar 
et al. 1997, 2000, and references therein), so it is worth revisiting the modeling to see if the de- 
viations of the average Wills et al. Fen models from the IZwl iron spectrum are now removed. 
Specifically, Wills et al. predict stronger Fen emission between ~1800 A and ~2300A, and weaker 
Fen right around Mgn for 3C273. In addition, the relative strength of various multiplets appear to 
differ between their models and the I Zw 1 spectrum. Variations in physical conditions between the 
NLSls, of which IZwl is considered the prototype, and the Wills et al. quasars may also explain 
some of the differences. 

Verner et al. (1999) present and discuss their numerical simulations of Fen emission spectra 
based on updated iron atomic data, iron line lists (many of which are also used here for line iden- 
tification), and photoionization modeling. They display the Fell emission spectrum for different 
densities, photon flux, micro-turbulent velocities below 100 km s -1 , and iron abundances. Although 
these line widths are much narrower than those of AGNs and quasars, complicating a comparison, a 
few differences between the synthetic spectra of Verner et al. and Wills et al. and our empirical Fe n 
and Fein templates can immediately be pointed out. (1) The synthetic iron models contain Fell 
transitions only, whereas the I Zw 1 spectrum also contain Fe III transitions. (2) The synthetic spec- 
tra (Wills et al. 1985; Verner et al. 1999) show a decrease in Fen emission above ~2700-2800A. 
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We are not able to see this in our template as contamination by the Balmer continuum (> A 2800 
A especially) is present. The HST spectrum of IZwl alone does not allow the Balmer contin- 
uum contribution to be sufficiently constrained. Complete photoionization modeling is required to 
properly account for it. However, the Balmer continuum contribution decreases continuously from 
~3600A down to ~2500A. Below ~2800A it is relatively weak, and it so should not significantly 
contaminate our 3100A-limited iron templates. 

Detailed comparison between the Fen template and the most recent Fen models by Verner 
et al. (1999) holds great potential for significantly improving our understanding of the iron emitting 
mechanism and mapping the physical conditions under which the iron emission is radiated. Hence, 
empirical iron templates remain useful even with the emergence of synthetic iron emission spectra. 

4. Application of the Iron Template 

4.1. Broadening of the Iron Template 

In order to fit and subtract iron emission from any quasar spectrum we need to broaden the 
iron template to match the line width of that spectrum. This is done by convolving the original 
template with Gaussian functions of different widths thereby creating a grid of spectra with a range 
of (quantized) widths. 

Three steps are necessary in the iron template broadening-process: (1) The standard deviation, 
Oconv, of the convolving Gaussian profile was first estimated using 

Vconv = FWHM cont ,/2v^2 = ^FWHMqso 2 - FWHMi Zwl 2 /2V2^2, (1) 

where FWHMiz w i = 900 km s -1 is the line width of the IZwl spectrum. (2) To apply a con- 
stant velocity broadening the computations were carried out in logarithmic wavelength space since 
d(log\) = dX/X = dv/c. The process of rebinning the broadened template back to linear units 
results in a small additional broadening. To ensure that the final broadened template has the 
desired width, we performed broadening simulations on artificial data consisting of a single Gaus- 
sian line feature of width 900 km s" 1 , where the width of the convolving Gaussian profile, a, was 
adjusted until the FWHM of the broadened feature, measured in linear wavelength space, matched 
the desired width to within our measurement errors (~lkms _1 ). (3) The iron template was then 
convolved (in logarithmic wavelength space) with a Gaussian profile with the adjusted value of a 
and rebinned back to linear units. 

We emphasize that because quasar spectral lines are often affected by blending, it is important 
to measure the resulting convolved width using a single artificial spectral line as outlined above, 
as opposed to estimating the width on the broadened (blended) template itself. We find blending 
effects to overestimate the line widths by ~15-25%, based on simulated line blends. Isolated 
narrow line cores, like Civ and Mgn are, however, expected to be well-determined. 
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4.2. Cleaning Quasar Spectra for Iron Emission Application of the Method 

In order to apply this UV iron template to quasar spectra, both the line widths and the iron 
emission strength of the template must be matched to that of the quasar spectrum. We follow 
the basic method of estimating the iron emission strength and subtracting the iron emission used 
by BG92 on the optical (AA 4250- 7000A) iron emission of their sample quasars. This method 
consists of artificially broadening the iron template to a number of widths by convolution with 
a Gaussian profile (§ 4.1) thereby generating a two-dimensional grid of templates with different 
line widths. BG92 scaled this "two-dimensional template" by a number of arbitrary, but fixed, 
strengths, hence creating a three-dimensional grid of templates (with dimensions: wavelength, line 
width, and line strength). The three-dimensional iron template was subtracted from the target 
spectrum (expanded to a data cube) and the best residual (one-dimensional) spectrum was chosen 
by manual inspection. 

Similar to BG92, we created a three-dimensional iron template consisting of five different 
scalings of a two-dimensional template, containing a number of line widths (in the range 1000 km s _1 
— 15 000 km s -1 in steps of 250 km s -1 is usually sufficient) which are broadened versions of the 
original iron template (intrinsic width 900 km s -1 ). We also performed a manual inspection of the 
residual iron-subtracted spectra. Our method, however, differs from that of BG92 in the following 
ways. After fitting and subtracting a power-law continuum fit 17 to the iron-contaminated quasar 
spectrum, we iterate to determine the optimum parameters as follows: 

• Determine the iron template (primary) normalization for a given quasar spectrum interac- 
tively from its iron emission strength in pure iron emission windows, specified by the user 
(examples are listed in Table 5). There is a separate normalization for each (1-D) template 
spectrum with a given broadening. 

• Subtract a (three-dimensional) data cube consisting of five scaled, two-dimensional templates 
from the target spectrum. The scalings are (five) fractions of the computed normalization, 
and can be chosen arbitrarily. A useful range is 50% — 150% of the normalization factor 
determined above, though the scalings /100% are rarely needed except for cross-checks. 

• Compute the residual flux and \ 2 of all the fits in the iron emission windows. 

• Manually inspect the residual spectra and pick the best (one-dimensional) iron-subtracted 
spectrum. 

• Add the previously determined continuum fit for the quasar back into the iron-subtracted 
spectrum and refit the power-law continuum. (In iron-contaminated spectra only a few nar- 



Based on our experience, optimal fitting of the iron emission in any quasar spectrum requires a power-law 
continuum fit, since a power-law continuum was initially subtracted from the I Zw 1 spectrum before generating the 
iron template. 
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row, pure continuum windows may exist. After a preliminary subtraction of the iron emission, 
larger continuum windows can generally be used). 

• Subtract this (new) continuum fit from the original target spectrum. 

• Determine the new normalizations as in the first item above and repeat all steps iterating 
over both continuum setting and iron emission strength until both fits (continuum and iron 
emission) converge and the final iron-subtracted spectrum is satisfactory. 

By estimating the template scaling factor interactively, we are not limited to an a priori set of 
quantized scalings (the quasar spectra can have any normalization and iron emission strength), 
and we ensure a good starting point for the iron spectrum fitting. Computations of both the \ 2 
and the residual flux in the user-defined iron emission regions help us to determine the best-fitting 
iron spectrum objectively. Iterating over both the iron emission strength and the continuum setting 
improves the fit to them both, especially in heavily contaminated quasar spectra where it may be 
hard to define pure continuum regions over which the underlying global continuum can be fitted. 

4.3. Sample Iron Emission Subtraction from Quasar Spectra 

We fitted the iron emission in four high-redshift and one low-redshift quasars, demonstrating 
the successful application of the UV iron template. We briefly summarize our main conclusions 
and then comment on these results for each individual quasar. The objects and their spectra 
were chosen to be representative of typically available data for the high-redshift quasars and to 
contain a range in Fe n emission strength. The 3C273 HST spectrum has the advantage of covering 
the entire UV range. The iron template is also successfully applied to the Large Bright Quasar 
Survey (Forster et al. 2001). We note that a perfect fit of the IZwl iron emission templates to 
the iron emission in other AGNs is never expected. This is owing to the strong dependence of the 
iron emission spectrum on the physical conditions in the emitting gas and on the BLR geometry 
(Netzer 1980). This is manifested in AGN spectra by variations in the strengths of the different 
multiplets among individual AGN (see discussion of 3C273 below). However, we find the template 
fitting to work well enough to be a valuable tool for eliminating and studying the iron emission in 
AGNs at least until accurate theoretical iron emission models are developed. 

A reasonable match to the iron emission in the quasar spectra presented here could be obtained 
by scaling and broadening the combined Fe n and Fe in UV templates, as described in § 4.2. Similar 
to the template, Fen emission is rather common around Mgn and between Civ and Cm]. The 
need to include one or more of the Fen UV191, Fein UV34, and Fein UV47 multiplets, which are 
strong in I Zw 1, in the applied template depends on the individual object. Improved iron template 
fits were generally obtained by sub-dividing the template and scaling each sub-spectrum (and hence 
selected groups of multiplets) separately. Good fits were often separately obtained for the Ly a — 
Civ, Civ- Cm], Cm] 2300A, and ~2300 3100A regions. 
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The most frequent and significant differences in strengths and multiplet ratios of the Fe II and 
Fe in transitions among the objects considered here are in the regions shortward of C iv, between 
~1900A and ~2300A, and in the Fell-bump around 2500A. The (undivided) template generally, 
but not always, overpredicts the A1400 - 1530 A iron emission in these quasars 18 . A separate 
(re-)scaling of this iron emission is acceptable because we have no empirical or simple theoretical 
constraints on the strength of this emission based on the iron emission, say, around Mgn. The 
variation in the A1900 - 2300 A emission strength is possibly due to differing amounts of dust 
extinction along the line of sight towards individual objects. The general weakness of the A2000- 
2300A emission in the template suggests that I Zw 1 may be subject to high levels of dust extinction 
perhaps from its host galaxy and consistent with the observed strong infrared emission (Rieke & 
Low 1972; Rieke 1978). The iron multiplet ratios may also (or instead) vary significantly in this 
region among quasars. This discrepancy is not unique to the quasars presented here (Figs. 6 and 7), 
but is also seen in other broad-lined quasars (e. g., Wills et al. 1985; Steidel & Sargent 1991; Corbin 
& Boroson 1996; M. Vestergaard et al. 2001, in preparation). Residual emission is similarly present 
in the sample iron-subtracted spectrum of Corbin & Boroson (1996; their Figure 2). 

The overall and important conclusion is that a subdivision of the ~2000 A wide Fe II and 
Fein UV templates may often be necessary because the individual multiplet strengths across the 
spectrum vary between objects. Similarly, the Fen and Fein templates may require different 
scalings as evidenced by the absence of Fein emission in some sources (e. g., Q1451+1017; Baldwin 
et al. 1996). Also, there is no strong direct coupling expected between the UV Fein and Fen 
emission. The Fen is emitted from the partially ionized hydrogen zone while the Fein originates 
in the fully ionized zone. Some correlation is expected to exist as both ions are (at least partly) 
excited by the UV continuum, but the current theoretical models are not accurate enough to firmly 
establish a simple scaling relation between the strong multiplets of UV Fen and Fein, if one exists 
(A. Pradhan 2000, private communication). Therefore, a separate scaling of the UV Fen and Fein 
emission templates is allowed in order to optimize the fit to the AGN iron emission. However, this 
was not necessary for the AGNs presented here. 

We first present the fitting results on 3C273 as its AA940 - 3200 A HST spectrum allows the 
full range of the template to be tested on a single AGN. 

Q1226+023; 3C273 (z = 0.157): The 3C273 spectrum (Fig. 6) does not contain any obvious 
features of the Fen UV191, Fein UV34, and Fein UV47 multiplets so the spectrum was fit with 
templates excluding these features. The combined Fen and Fein template did a good job, once 
subdivided into five independently scaled segments (at AA 1540, 1912, 2423, and 2656A). A scaling 
of the full range template to the iron emission redward of C m] and in the small blue bump yielded 
consistent flux levels blueward of Civ if no Niv] A1486 line is present. We could not distinguish 



18 This may partly be due to the fact that I Zw 1 brightened (§ 2.1) between the G130H and the G190H and G270H 
spectra were observed, and that the iron spectrum does not simply scale with the continuum level. 
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between this fit and one with a separate scaling and a Niv] line (the latter fit is shown in Fig. 6). 
The fitted iron model in the range between C iv and C in] yields residuals coinciding with the 
expected wavelengths of known lines, such as Hen A1640, Olil] A1664, Alii A1670, Nlli] A1750, 
Sin A1814, and Aim A1857 yielding confidence in the template and the model fit. We found the 
Mgn profile in the iron subtracted 3C273 spectrum to show a slightly blue asymmetric profile. This 
profile could be fitted very well with two Gaussian functions with the same FWHM of 3800 km s _1 
(one shifted — 3400 km s -1 from the narrow peak of Mgn; FWHM is consistent with the widths of 
the UV lines and the narrow Balmer lines, ~4000kms~ 1 ). This gives us confidence that the iron 
template in the Mgn region is representative of the average Fen bump in AGNs. 

The 3C273 spectrum shows stronger iron emission than I Zw 1 between the C iv and C III] lines 
and in the AA2000 — 2300 range (Fig. 6) relative to the remaining UV iron emission. As the I Zw 1 
template does not contain much iron emission in the AA2100 - 2250A range, no further attempt to 
model this additional emission is possible with the current template. The I Zw 1 template clearly 
overestimates the 3C273 iron strength around 2500A. The emission is rather flat across this region 
contrary to that in the template. Two different attempts to model this region were made (Figs. 6b 
and 6c). One 'model' fits the blue part while the other fits the red part, but neither does a good 
overall job. The original template has much stronger emission at 2500 A than that shown in Fig. 6. 
In the case of the ~2100A and ~2500A 'bumps' a different template is clearly needed to account 
for this emission. Redward of Mgn the template shows a similar shape as the observed Fen but 
slightly underestimates the emission. It is likely due to a slight overestimation of the underlying 
continuum level in I Zw 1 owing to the truncation of its spectrum at 3089A rest frame by the FOS 
G270H grating. 

The power-law 19 continuum has F>,(1500 A) = 1.608xl0 -13 ergs cm~ 2 s _1 A -1 and slope, 
a\ = 1.67 (a u = 0.33). The original spectrum is shifted by +1.5x 10~ 13 ergs cm -2 s _1 A -1 in Fig. 6a 
and by +0.25 x 10 -13 ergscm~ 2 s _1 A -1 in Fig. 6b. In this figure the residuals from the 'blue 2500A' 
iron model, shown in the middle, represents the original flux level of the spectrum. The residuals 
from the 'red 2500A' iron model, shown as the lowest of the three spectra in panel (b), is shifted 
by —0.25x10 ergs cm" gure 6c shows the two 2500A iron models more clearly; the 

spectrum and the models are continuum subtracted here. The iron model has a FWHM of 4000 
kms -1 . The processing of the HST data is outlined in Appendix B. 

Q0020+022 (z = 1.798): There is no clear indication of strong iron emission in the A1250- 
1500 A range and the continuum level is relatively lower than the flux level in the A1600-1900A 
range. The flatness of the A1930-2050A region also indicates very weak iron emission, so this 
region was adopted as a continuum region. The power-law continuum has normalization, F>,(1500 
A) = 1.753xl0~ 16 ergs cur 2 s -1 A -1 and slope, a A = 1.94 (^ = 0.06). The combined Fen and 
Fein template was fitted to the data between A1530 and 1900A only (Fig. 7a) yielding residuals 



F v ~ v- a " and F A ~ \~ a *;a x = 2 - a„ (e. g., Weedman 1986). 
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which coincide very well with the expected positions of Hen A1640, Oiii] A1664, Alii A1670, Niii] 
A1750, and Sin A1814A. The residual spectrum is shifted by — 0.5xl0~ 16 ergs cm -2 s^ 1 A~ x . The 
iron model has FWHM of 6500 km s -1 . 

Q0252+016 (z = 2.457): An improved fit to the iron emission (Fig. 7b) could be obtained by 
including the FenUV191 A1786 and FemUV34 AA1895, 1914, 1926 features in the templates and 
by separately fitting the A1400 - 1550A emission which is weaker than in the UV template, but still 
significant. It is unclear whether weak Fen emission still remains in the blue wing of Civ A1549, 
though both wings of this line indicate the presence of a somewhat broad underlying component. 
The power-law continuum has Fa(1500 A) = 2.807xl0 -16 ergs cm -2 s _1 A -1 and slope, a^ = 1.28 
(oj, = 0.72). The residual spectrum is shifted by — 1.0xl0~ 16 ergscm" "2 s- 1 A" 1 . The iron model 
has FWHM of 5000 km s" 1 . 

Q1629+680 (z = 2.478): The iron emission is relatively weak (Fig. 7c) compared to the other 
quasars presented here, as indicated by the lower emission level between Civ and Cm]. Based on 
the C iv line width of FWHM~4000 kms -1 , a reasonable but weak model fit could be made using 
the combined Fen and Fein templates in which Fen UV 191 A 1786 and Fein UV 34 AA1895, 1914 
were excluded. One exception is the region between A1400 and C IV, where the UV template strongly 
overpredicts the Fen emission. In fact, the data does not indicate noticable Fen emission there. 
Thus, the iron emission was fitted using a version of the template in which the A1400-1550 A iron 
emission is excluded. The power-law continuum has Fa (1500 A) = 8.903 xl0~ 17 ergs cm -2 s -1 A -1 
and slope, ot\ = 2.3 (a v = —0.3). The original spectrum is shifted by +1.0 xlO -16 ergscm -2 s -1 A -1 . 
The iron model has FWHM of 4000 km s _1 . 

Q2345+061 (z = 1.540): A scaling and broadening of the combined Fen and Fein templates 
(with Fell UV 191 A 1786 excluded) to ~2500kms~ 1 offer a reasonable overall fit to the iron line 
emission (Fig. 7d). Hence, no sub-division of the UV template was performed. The general shape of 
the iron emission longwards of Mg 11 is reproduced by the template but a slight drop in the emission 
level is seen in the quasar data owing to the enhanced calibration uncertainty in this region where 
second order light merges with the primary order spectrum. The main bumps and wiggles in the 
Fe 11 bump around Mg II can be reproduced with the exception of the A~2300 - 2500 A region where 
excess emission is clearly present. Including the Fein UV47 AA2418, 2438 feature in the template 
(shown in Fig. 7d) does not fully account for the iron emission in this region. Excess emission 
is also clearly present at A~2050-2150 A. A good coincidence of the iron strength is found for 
~1800-2000 A, based on a scaling of the template Fen bump around Mgn to the data, inspite of 
the larger noise level below 2100 A. The large noise levels at A< 1800 A renders the discrepancies of 
the data with the template fit irrelevant. The power-law continuum has Fa(1600 A) = 3.62xl0~ 16 
ergscm -2 s" 1 A -1 and slope, ot\ = 1.72 (a v = 0.28). The original spectrum is shifted by +1.0xl0~ 16 
ergs cm ^s -1 ! -1 . The iron model has FWHM of 2500 kms . 
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5. Discussion of the I Zw 1 Spectrum 

The spectrum of I Zw 1 is rich in emission lines, and in particular this high S/N HST spectrum 
reveals a plethora of weaker lines (L97; Fig. 1) not usually seen in quasar spectra. The spectrum has 
strong low-ionization lines while some of the higher-ionization lines are weaker than expected for 
an average quasar (Francis et al. 1991; Zheng et al. 1997; Wilkes 2000; L97 and references therein; 
see discussion on Siiv+Oiv] and Civ below). Both [Om]A5007 (BG92) and Cm] A1909 (this 
spectrum) are weak. 

While a detailed emission line study is not the subject of this paper, we briefly discuss our 
results in relation to earlier work on this object. 

5.1. Line Parameters 

5.1.1. Line Widths 

Our measurements and line fitting show that the widths of the fitted iron emission components 
are rather narrow; ranging from ~300 to 900 kms -1 . E. g., the narrow line core of FellUV191 
A 1786 has FWHM of 550 km s" 1 and a similar fit to Fe in UV47 A 2418 yields a FWHM of 715 km s" 1 
(Table 4). This narrow width may be an intrinsic property of the iron emission mechanism or an 
effect of the continuum level uncertainty (the line width will appear narrower if the continuum is 
set too high). Our narrowest line measurements of ~300-400 kms -1 of the weaker non-iron lines 
are consistent with the findings of van Groningen (1993), who concluded, based on work by himself 
and Phillips (1976), that IZwl has three emitting regions, one of which has a velocity dispersion 
of -400 kms" 1 . 



5.1.2. Line Velocity Shifts 

The peaks of the broad emission lines are blueshifted (Table 3) relative to the systemic redshift 
(z em = 0.061), defined by the Mgn doublet (Appendix A) and Hi radio measurements (Condon 
et al. 1985). The high ionization lines have a higher blueshift (average shift ~1540±500 kms -1 ) 
than the low ionization lines (~500±270kms~ 1 ; see also Table 3, L97, and below for exceptions). 

A significant amount of blueshifted emission is evident in lines such as N v, the C iv+Si n blend, 
and He II A 1640 (Figs. 5a, 5c, and 1, respectively). A special case is that of Hen where practically 
all the emission is blueshifted, reaching ~1300— 2200 kms -1 relative to the quasar restframe (see 
also L97). The Nv line profile can be fitted as two emission components originating at different 
velocity shifts, ~ -900 kms" 1 and ~ -2400 kms" 1 (Fig. 5a, Table 3). The ~900kms~ 1 blueshift 
of Nv is similar to that of the narrow Civ component, the broad Aim component, and one of 
the narrow components of Ly a and Siiv+O iv], respectively. Most of our line shifts are consistent 
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with those of L97 to within the errors and the spectral resolution (Leitherer 1995; Tables 3 and 4). 
We suspect the velocity shift differences < 1300 km s -1 with L97 for Sin A 1263, Sin+0 1 A 1306, 
SilV+O iv] A 1400, and Niii]A1750 can be explained by a combination of uncertainties in the 
line positions of < 0.5— lA (i. e., 250— 500 kms -1 ; conservative error), resulting from the heavy line 
blending present in both studies, and the slightly degraded resolution in our G130H+G190H specta 
(§ 2). 

L97 suggested an outflowing component is responsible both for the blueshifted line peaks (of 
both low and high ionization), the blueshifted (~2000kms~ 1 ) emission (blue wing asymmetry) and 
the associated weak UV absorption in Lya, Nv, and Civ which they detect. 



5.1.3. Line Strengths 

L97 use the doublet ratios of Mg n and Al in to deduce the location of their emitting regions 
relative to the 'outer BLR boundary'. As discrepancies are seen between our line measurements 
and those of L97, it is of interest to briefly discuss the reasons and the implications for their BLR 
size results. At A < 1500 A the main differences are due to the fact that L97 does not scale the 
G130H spectrum to match the level of the G190H and G270H spectra, as done here. 

The line measurements at A > 1500A are different in part due to different Gaussian fitting 
techniques (cf. Laor et al. 1994), but mostly due to L97's use of a local continuum (A. Laor, 
1997, private communication) as opposed to a global one, as done here except at the Fen bump, 
A2300 - 3090 A. The EW measurements in regions of isolated and unblended emission lines (where 
the continuum is well determined) agree well (CnA1335, A 1345-feature, and Sin A 1814). How- 
ever, for the remaining lines (see Table 1 in L97 and Table 3, this work), the EW measurements 
differ by —75% to 180% and the line fluxes by —70% to 215%, where the faint or blended lines 
deviate the most (|AEW| > 30% and |A(line flux)| > 30% in lines such as Sin A 1260 [blended], 
O ill] A 1664, Al II A 1670, Fen A1786, Aim AA 1854, 1863, Nil] A 2141, Cn]A2326). These discrep- 
ancies measured from the same data illustrate the significant uncertainties associated with line and 
continuum fitting in quasar spectra. 

We can not easily confirm the suggestion by L97 that Al in AA 1854, 1863 can probe the BLR 
size, as it appeared only marginally thermalized (ratio = 1.25:1). Our measurements of the narrow 
Aim doublet 20 (ratio = 0.9:1) suggests it is thermalized, possibly being emitted somewhat closer 
to the continuum source, i. e., inside the BLR outer radius. We also find Mgn to be (entirely) 
thermalized (doublet ratio = 1:1); L97 measure a ratio = 1.2:1. Contrary to L97 we find a much 
better overall fit to the narrow Mg II emission using more than two components; four components 
are required to ensure an equal line width in each doublet (Table 3). Thus it appears that neither 
Mgn nor Aim is a suitable probe of the outer boundary of the BLR (see discussion by L97). 



The Al in doublet components have equal widths to within the spectral resolution. 
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5.2. The C in] Complex 

The density-sensitive lines, C ill] A 1909 and Si m] A 1892 are important BLR density diagnos- 
tics. At first sight their relative strengths in IZwl are highly unusual: Cm] is ~10 times weaker 
and Sim] significantly stronger 21 than usual leading L97 to argue for unusual high densities (~10 11 
cm -3 ). However, in reality this comparison is difficult to make. Published measurements of Cm] 
line strengths (e. g., Wilkes 1986, Francis et al. 1991, Baldwin et al. 1995) are generally measure- 
ments of the entire Cm] complex including Sim], Aim, and Fein emission, which are all strong in 
IZwl. In typical quasars with much larger intrinsic line widths these lines are so heavily blended 
that detailed deblending techniques often cannot uniquely determine the relative contributions from 
individual transitions. Note, if Si m] and Fe m are both strong at high densities the complex would 
not necessarily appear asymmetric at FWHM > 3500 km s -1 . By artificially broadening the IZwl 
spectrum to FWHM of 3000 km s -1 and 5000 km s -1 , a typical range for quasars, we estimate the 
C m] line strength measurements (line flux and EW) in a similar, broad-lined source would be 
overestimated 22 by a factor of ~2. 

Detailed modeling of the full UV spectra of other NLSls also concludes that these sources have 
high density emitting regions (~ 10 11 — 10 12 cm~ 3 , Kuraszkiewicz et al. 2000), lending support to 
L97's conclusions. Given that the strong Fe m UV34 multiplet is severely blended with the C in] and 
Sim] emission lines, we note that it is likely that both L97 and Kuraszkiewicz et al. overestimate the 
strength of the Si m] line as they underestimate the strength of the Fe m UV34 1895A transition. 
The combination of a resulting smaller Si III] / C III] ratio and strong Fe III and Al III contributions 
will in the framework of photoionization models (e. g., Fig. 3e, Korista et al. 1997) then further 
strengthen their conclusions that the BLR densities, at least in NLSls, are very high. 

A theoretical estimate of the relative triplet strength for the Fein UV34 AA1895 1914, 1926 
transitions for AGN physical conditions is not readily available. The relative significances of the 
possible excitation mechanisms are not fully established in part due to the complexity of the 
Fen and Fein emission. The Fein UV34 excitation mechanism is probably not due to electron 
impacts owing to the high energy levels of the multiplet transitions (several eV) compared to the 
(typical) plasma temperature of ~10 4 K (~ 1 eV). The multiplet could be due to fluorescence and/or 
photoionization/recombination, as is likely to be the case for the UV Fein emission in general (A. 
Pradhan 2000, private communication). In the optically thick regime all the transitions should 
have equal strengths due to thermalization (e. g., L97). Hartig & Baldwin (1986; hereafter HB86) 
estimate the relative Fe m UV 34 multiplet strengths by fitting the AA1895, 1914, 1926 transitions in 
the spectrum of H0335— 336 where this multiplet dominates the Si m] and C m] emission. Assuming 
H0335— 336 has no Sim] and Cm] emission at all, they find a multiplet ratio of 0.9:1:0.3, while 
including Sim] and Cm] in their fits yields a relative flux ratio of 0.9:1:0.7. The IZwl spectrum 



The critical density of Si ill] is 1.1 xlO 11 cirT 3 (Baldwin et al. 1996; L97). 
The Aim line is excluded in the comparison. 
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is consistent with either of these ratios, as we discuss next. 

In order to estimate the likely contribution to the Cm] complex from the Fein UV34 triplet, 
we modeled this multiplet with special emphasis on the relative transition strength, on varying its 
contribution and noting the effect on the Si III] and C in] model fits (see Table 6) . We adopted the 
approach by HB86, who use the emission profile of Fen UV191 to model the Fein triplet (J. A. 
Baldwin 2000, private communication). We note that due to the strong blending in the Cm] line 
complex and the non-orthogonality of Gaussian functions, the model fits are not unique. This is 
also clear from the fact that a range of reasonable model fits to the Fein UV34 emission can be 
made. Assuming there is no Si III] and C m] emission at all (or broad emission thereof at least) a 
triplet ratio of 0.91:1:0.46 can be fit (not shown). However, the significant residuals at A 1890 and 
A 1907 strongly argue that Sim] and Cm] emission is present. Several fits were made both where 
the triplet ratio was fixed at 1:1:1 and where it was allowed to vary freely. Sample fits are shown 
in Figure 5e and Table 6 (Note that the triplet is located at AA1893 1912, 1924). The best fits 
were found for triplet models where the A 1914 A transition is the strongest. We cannot distinguish 
between the various fits nor between the optically thick and thin cases. More advanced modeling 
is required. A reasonable approach in fitting the Fein in this region until more is known about 
the UV34 triplet may be as follows. If Aim is strong, Fein (including UV34) is highly likely to be 
present and similarly strong (HB86). If so, two fits can be made: one ignoring the presence of Fe m 
UV34, and one in which its contribution is maximized. This allows an upper limit to the fitting 
uncertainty involved to be estimated. 

In Figs. 2, 3, and 5d the Fein template is shown with the UV34 triplet ratio 0.375:1:0.425 
(model B) for illustrative purposes only. We note that other strong Fe III transitions may be present 
in the C m] complex in addition to Fe ill UV34, as indicated in Figure 4. 

5.3. The Si iv+O iv] A 1400 and C iv Features 

The Siiv+Oiv] A 1400 feature is strong relative to Civ (A1400/Civ ~1.3; Table 3), and the 
Civ/Lya ratio (0.11) is also low compared to the typical values for quasars: 0.1-0.6 and 0.2- 
0.6, respectively (Wilkes 1986; Francis et al. 1991; Baldwin et al. 1995). We note that A1400/Civ 
measured from the non-scaled G130H spectrum (§ 2.1) is also unusually strong (~1). The line lists 
by Nahar (1995) and Nahar & Pradhan (1996) indicate a (conservative) upper limit to the iron 
emission contribution to the A1400 feature of ~20%. Thus, the A1400/Civ line ratio indicates a 
density n e > 10 n cm~ 3 according to the models by Rees, Netzer & Ferland (1989), consistent with 
the estimate by L97 of n e ~ 10 11 cm~ 3 and typical NLS1 densities (Kuraszkiewicz et al. 2000). At 
high densities a simultaneous strengthening of the Silicon lines and weakening of the Carbon lines 
occurs (e. g., Rees et al. 1989). 

Another contributor to the large A1400/C IV intensity ratio may be an ionization effect due to 
the unusually red spectral energy distribution of IZwl [e. g., similar to the explanation by Zheng 
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& Malkan (1993) of the Baldwin effect (Baldwin 1977)]. The luminosity brightening relative to 
the IUE data measurements (§ 2) appears to have increased these line ratios in agreement with 
Zheng & Malkan (1993), with the C iv complex flux increasing the least (~75 %; Ly a increased by 
~110%; Wu et al. 1983; Wang et al. 1996; Rodriguez-Pascual et al. 1997). 

The broad component of Civ is relatively weak compared to the other UV lines (Table 3; 
Fig. 5c). Marziani et al. (1996) define a strong broad C iv component in the i/STspectropolarimetry 
data 23 of IZwl, but they fitted the component to emission which we identify as Sin, possibly 
Fe II and blueshifted C iv. Given the lack of deblending or detailed line identification their broad 
component fit has an EW = 21 A, compared to ours of ~9 A. 

5.4. Is Iron Emission Associated with Outflows? 

Baldwin et al. (1996) study the nature of the kinematic components in the emission line spectra 
of seven quasars and, based on one of them, Q0207— 398, propose that three different components 
are present in the line profile. They argue that Q0207— 398 is a "misaligned" broad absorption line 
(BAL) quasar as the blueshifted line emission in Q0207— 398 is typically absorbed in BAL quasar 
spectra. Baldwin et al. also connect this blueshifted emission and absorption to the expanded 
photospheres of (bloated) stars close to the central source and argue that Aim A 1857 emission is 
an indicator of the presence of these stars. 

The UV spectrum of Q0207— 398 bears a remarkable resemblance to that of IZwl (narrow 
line cores, strong Aim, Fen, and Fein emission, high density emitting regions, and blueshifted 
emission in the high ionization lines), and the different kinematical components in Q0207— 398 are 
similar to the emitting regions at different redshifts in IZwl (Phillips 1976; van Groningen 1993), 
providing further support for the argument that IZwl is (also) a "misaligned" BAL quasar (L97). 
Baldwin et al. (1996) connect the Aim emission to outflows (and stars). So as the Aim and Fein 
emission strengths appear to be connected (e. g., HB86) we speculate that the iron emission itself 
is somehow related to the presence of the outflows and/or to the stars. The fact that the iron 
line profiles do not appear asymmetric or are significantly blueshifted indicates that this emission 
does not originate in the outflowing gas itself. In any case, strong iron (Fen and Fein) emission 
may be connected with high densities (HB86; Joly 1991; Baldwin et al. 1996; Lawrence et al. 1997; 
Kuraszkiewicz et al. 2000). 



Note, their flux level is different from our spectrum and that of L97. 
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5.5. Is the Iron Spectrum of I Zw 1 Typical? 

An underlying assumption for using the I Zw 1 spectrum as an iron template is that the iron 
emitting mechanism in this target is similar to that of typical quasars and that the spectrum 
provides a good representation of the iron transitions and iron line ratios observed in AGN and 
quasars. We therefore address the question of how typical the emission lines, including iron, in 
I Zw 1 may be. 

I Zw 1 has a number of unusual properties in addition to strong Fe n emission. It has strong Ca II 
emission at AA 8498, 8542, 8662 (the infrared triplet, permitted lines) and AA 7291, 7312 (forbidden 
lines) (Phillips 1976, van Groningen 1993, Persson & McGregor 1985), strong Cm* A 1175 emission 
(Laor et al. 1997a; L97), strong Si III] A1892, weak Cm] A 1909, strong IR emission (e. g., Rieke & 
Low 1972; Rieke 1978) plus strong Fein emission transitions in its spectrum (Figs. 1 — 3). 

Strong Can emission is observed in ~30% of all quasars and AGN (Netzer 1990) and is 
generally thought to be emitted in the deepest interior of BLR clouds (Persson & McGregor 1985; 
van Groningen 1993) with very high column densitites (N# > 10 24 - 5 cm~ 2 ). Strong Can emitters, 
however, do not otherwise appear different from those of the average AGN population (Netzer 
1990). 

Cm*A1175 is not a commonly detected metastable transition of Cm in AGN and has previ- 
ously most often been detected in absorption (e. g., Bromage et al. 1985, Kriss et al. 1992), perhaps 
due to the difficulty in identifying very broad, but weak features (§ 1.3). See e. g., Laor et al. (1995) 
and Hamann et al. (1997) for weak and marginal detections in emission. The line may result from 
the enhanced density (see discussion by L97) of the emitting medium already deduced for IZwl. 
Several studies show a trend toward strong Fe II, Fe III, and/or Al III at higher BLR densities (§ 5.4). 
If high densities are a common property of strong iron emitters, the iron spectrum is not expected 
to be unusual. 

Strong Fein features are more prevalent in IZwl (Table 4) than in previous AGN studies, 
although Fein features at 2070 A and ~2420A have been reported in the past (HB86; L97). Our 
study suggests, however, this difference is not real but due to misidentifications in the past. Francis 
et al. (1991) note unidentified features in the AA2000-2200A and AA2900-3200A regions, consistent 
with Fe III lines according to our identifications (see references to Table 4). The presence of both ions 
is potentially important for deriving the physical conditions in the iron-emitting regions, because 
the details of the iron emission spectrum are highly sensitive to the physical conditions (e. g., Netzer 
1980; D. Verner, 1997, private communication). 

Lipari et al. (1993) and Lipari (1994) link extremely strong optical Fen emission to starburst 
activity. Though I Zw 1 is currently undergoing vigorous star formation (based on Barvainis et al. 
1989; Eckart et al. 1994; Sanders & Mirabel 1996), the iron emission observed in AGN is not 
characteristic of emission from star forming regions. Pure starburst galaxies (e. g., NGC 7714) do 
not emit permitted Fell emission, only IR forbidden lines (e. g., [Fell] 1.6 ^m; L. Ho, 1998, private 
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communication) mainly because the stellar ionizing continuum is too soft. To obtain the permitted 
lines the hard, non-thermal continuum from a central AGN source is needed to penetrate to the 
high-density interior of the BLR clouds (e. g., Netzer 1990). 

Based on the above discussion, we find no strong evidence that the iron emission in I Zw 1 is 
unusual compared with that of the general quasar population. This is confirmed by our successful 
fitting of the iron emission in several quasar spectra (§ 4.3). 

6. Summary and Conclusions 

We have presented a UV iron template based on HST archival data of IZwl, and described 
the method with which the template was generated. Compared to previous empirical templates 
(Corbin & Boroson 1996) this covers a large range (AA1250 - 3090A), was generated by careful 
fitting of the non-iron emission and absorption lines, and allowed the generation of separate Fell 
and Fe in templates. We have demonstrated its application to fit and remove the Fe II (and Fe III) 
emission in spectra of several quasars, including 3C273, allowing subsequent studies of weak and 
heavily blended emission lines (e. g., Cm]) free from the large uncertainties otherwise associated. 
This shows that the iron emission in I Zw 1 is sufficiently similar to that in other broad-lined quasars 
to be useful. Although it has limitations, the iron template fitting process is an important tool to 
eliminate and study the iron emission in active galaxies, at least until accurate theoretical models 
are developed. 

We confirm previous results reporting the presence of blueshifted emission (i. e., blue profile 
asymmetry) in the spectrum of IZwl, especially in the high- ionization lines (L97), and of several 
emission regions of different (absolute) redshift, line width and physical conditions (Phillips 1976; 
van Groningen 1993). The presence of these regions along with the UV spectral properties arc 
consistent with the picture suggested by Baldwin et al. (1996) in which the blueshifted emission 
originates in the same region giving rise to the BAL troughs in BAL quasars. This may connect the 
NLSls with BALs, perhaps through source orientation. We also argue that (strong) iron emission 
may be connected with high densities and to outflowing material. 

We discuss the fact that the Siiv+Oiv] A 1400 feature is very strong relative to C iv. This can 
be explained by a simultaneous weakening of the Carbon lines and strengthening of the Silicon lines, 
common to Narrow Line Seyfert Is (Kuraszkiewicz et al. 2000), probably due to a combination of 
high densities and a low ionization parameter. 

Based on earlier IUE measurements we find I Zw 1 to have brightened a factor ~2 in its con- 
tinuum emission since 1978-1982. 
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A. Comments on Individual Spectral Features in I Zw 1 

Unless otherwise noted, the identifications of Fe n and Fe in emission features are made using 
the iron line lists available on the web and from recent publications (Fuhr et al. 1988; Penston et al. 
1983 ; Ekberg 1993; Giridhar & Arellano Ferro 1995 ; Nahar 1995; Nahar & Pradhan 1996; Quinet 
1996; Quinet et al. 1996; Kurucz & Bell 1995). 

When (iron) residuals are said to be removed or excluded (i. e., subtracted) from the template, 
it is done in order to prevent an overestimation of the iron emission strength and subsequent 
overcorrection for the iron emission when the template is applied to quasar spectra. In addition, 
some line fits may require a local continuum level (different than the global continuum setting). 
However, the EW measurements are always with respect to the global continuum. 

AA1120— 1135 emission: Fein UV1 emission is clearly detected, but is heavily affected by ab- 
sorption (Figs. 1 and 2). Due to the uncertainty in correcting for this absorption in the noisier end 
of the spectrum, we chose not to include the region A1075— 1135A in the template. 

C III* A 1176 emission: This emission line (Figs. 1, 2, and 5a) is an excited meta-stable level of C in 
and rarely detected in AGN spectra [but see Laor et al. (1995); Hamann et al. (1997), and see for 
absorption detections Kriss et al. (1992); Bromage et al. (1985)]. Its implications are discussed in 
§ 5.5, and Laor et al. (1997a, 1997b) also discuss possible excitation mechanisms of this transition. 

AA1150-1245 range: Faint emission of Fen is expected in this wavelength range (Figs. 1 and 2). 
The leftover residuals from subtracting the fits to the very strong, non-iron lines (especially Ly a ) 
are rather noisy and do not obviously match expected Fen transitions. We therefore choose to 
remove these residuals from the template (Fig. 3). 

Ly a A1216 emission: Two narrow components are necessary to reproduce the shape of the narrow 
line core (Fig. 5a). The weaker of the two components is blueshifted (~900kms~ 1 ) relative to the 
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peak position of the stronger component. The high S/N data constrain the fit well. Although the 
Lya profile does not closely resemble that of the Balmer lines (cf. Ha and H/3 modeled by BG92, 
and L97), the profile shapes of all three are well reproduced by one broad and two narrow Gaussian 
components, one of which is blueshifted resulting in an asymmetric profile. 

O vA1218 emission?: The feature clearly appears once the fit to Lya is subtracted and the line 
position is well defined at 1218A (Fig. 5a). The identification (Penston et al. 1983) is uncertain as 
one would expect a high blueshift of the line similar to the other high-ionization lines. 

N v AA1238, 1243 emission: When the absorption is corrected, a blue asymmetric profile (Fig. 5a) 
is left, similar in appearance to the non-deblended profiles of Siiv+Oiv] and Civ. The doublet 
profile is fitted well with three Gaussian components, one of which is considerably blueshifted 
(~2500 km s _1 ) relative to the A1243 line. The remaining doublet feature is blueshifted ~900 km s _1 
similar to the blueshifted Ly a emission, relative to the rest frame. 

AA1230, 1234 absorption: L97 identify this absorption (Fig. 5a) as due to (~2000kms~ 1 ) 
blueshifted Nv absorption intrinsic to IZwl, supported by their detection of weak associated 
absorption in Lya and Civ also at this blueshift. 

The absorption lines cannot be due to OiA1306 and Sin A1309, although the wavelengths coincide 
very well with their laboratory wavelengths, as these transitions originate in slightly excited states 
which are not common in the ISM (Savage et al. 1993). 

A1257 — 1268 emission: This feature is identified as SinA1263 (Fig. 5a). Due to the detection 
of Fen UV9 emission at AA1270-1280, Fen (UV9 AA1250-1280A; Fig. 4) is also expected to be 
blended with this Sin multiplet. The wavelengths of the Gaussian components (Table 3) fitted to 
the Si ii line blend are uncertain, due in part to the short wavelength emission being blended with 
the Lya fit (and hence partly subtracted) and in part to blending with Fen. 

Sim* A1297 emission?: The identification of this feature is not confirmed (e. g., Fig. 2). It 
is not clear whether the emission is due to Sim* emission or to a blueshifted component of 
OiA1306 + Sin A1309, similar to that found in the high-ionization lines. In the latter case the 
blueshift is ~2000kms~ 1 , in agreement with the blueshift of the high-ionization lines. A third 
possibility is faint Fe in A1298.6 and Fe n AA1297, 1299 emission which is expected based on the iron 
line lists, although the observed emission feature appears too strong. 

OiA1304 emission: Faint Fen appears, especially in the blue wing, but the iron emission could 
not be deblended (Fig. 2). 

A1343 feature: The identification is uncertain (Figs. 1 and 2). No obvious Fell, Fein or other 
commonly observed quasar broad emission line matches the wavelength position. Candidate identi- 
fications are Oiv, Can, Fein and Fev. A coincidence occurs with some Fein transitions, as noted 
in Table 3, but it is not clear whether an isolated feature, as observed, is expected. Can is also 
observed in the IR (§ 5.5), but this A1343 feature is redshifted with respect to the expected Call 
position (Table 3). 
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Siiv+Oiv] A1400 emission: Faint Fen emission is expected at A1393 and A1408. A faint Fein 
feature is expected at A1395, which is identified in the residuals after subtracting the Siiv+Oiv] 
fits. The strength of Fein A1395 is likely to be underestimated: the feature is heavily blended with 
the Siiv+Oiv] emission and the data do not permit sufficient constraints to be placed on each 
contributing emission component. 

Blueshifted emission is clearly present in this line complex (Fig. 5b). A deblending of each of 
the components is not straight forward due to the severe blending of the components emitted in 
the restframe with the blueshifted line emission. Hence, each individual component has a larger 
uncertainty than the complex fit as a whole. We interpret the complex as follows: the "restframe" 
emission has a general blueshift of ~900 km s _1 relative to the Mg n line peak (consistent with Ly a , 
Nv, Civ, and Aim) and the excess blueshifted emission has a projected velocity of ~ -1100 km s -1 
relative to the former ("restframe") emission and a total blueshifted velocity of ~2000kms~ 1 
relative to the true restframe. The complex clearly has a broad component, which displays the 
strongest blueshift of 2000 km s -1 (Table 3). 

The line complex as a whole is strong relative to Civ, most likely because Civ is weak and/or 
the density is high and the ionization parameter is low (see discussion in § 5.3). Alternatively, the 
iron emission may be stronger than we currently suspect, making the Siiv+Oiv] emission appear 
stronger and contributing to the high A1400/C iv line ratio. 

Niv] A1486 emission: This emission line is not detected (Figs. 2 and 5c). It may, however, be 
weak and heavily blended with iron emission. 

AA1515, 1517 dips: These features (Fig. 5c) appear very similar to some of the Galactic absorption 
doublets, but comparison with the low dispersion IUE spectrum (Lanzetta et al. 1993) clearly show 
that this dip is a lack of emission surrounded by line emission, which we identify as Fe in (Table 4). 

Si II AA1527, 1533 emission: This doublet transition is almost certainly present (Fig. 5c). The 
lines are, however, strongly blended with C iv and possibly iron emission, and the individual line 
profiles and their strengths cannot be determined with confidence. Due to the apparent weakness of 
the broad C iv component we suspect the contribution from the Sin doublet to be significant. An al- 
ternative interpretation is that C IV has a strong and rather blueshifted component of ~6000 kms" 1 
relative to the restframe. In that case, there must be several blueshifted emitting regions as the 
observed emission does not appear to reflect the doublet flux ratio which is estimated for the narrow 
Civ components (Table 3 and Fig. 5c). Fein UV84 and weaker Fein features are also expected to 
be present (see below). 

C iv AA1548, 1550 emission: Blueshifted emission appears present (Fig. 5c) but, due to the heavy 
blending, a well constrained Gaussian component fit (deblending) to determine the blueshift and 
flux is not possible. The data only support a weak broad component (Fig. 5c; Table 3). Narrow 
profiles were fitted to the peaks of the Sin AA1527, 1533 and Civ AA1548, 1550 doublet lines which 
are clearly resolved. The best fit to the Civ narrow emission has slightly different widths of the 
transitions (Table 3). If the line width is fixed and a flux ratio of 2:1 is approximated a Gaussian 
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residual appear (see the alternative fit to the narrow C iv emission in Table 3). In addition, a two- 
component feature at A1536 and A1539, coinciding with Fe in UV84 emission (see below and Fig. 4), 
is required to fit the complex well. An alternative fit has two narrow components with a flux-ratio 
of approximately 2:1 (Table 4) suggesting an additional blueshifted (~2300kms~ 1 ) component of 
C iv emission. The likelihood that this emission is Si I is probably rather low as expected strong 
transitions at ~A1542, A1547, A1553, and AA1550-1600A (Morton 1991) are not seen. 

Fein A1520 - 1540 emission: The UV84 multiplet of Fein among other Fein lines of detectable 
strengths (AA 1515.5, 1516.2, 1518.8, 1524.5, 1525.0, 1526.0, and 1527.0; see Table 4 and Fig. 4). are 
expected close to the positions of Sin and Civ; the most probable lines are expected at AA1526, 
1527, 1531-1532, 1538-1539, and 1550-1551 (Nahar & Pradhan 1996). If these transitions are 
present Sin and Civ may be even weaker relative to the other resonance UV lines. The IUE 
spectrum presented by Lanzetta et al. (1993) does not show any signs of strong Fein emission in 
this region. But since IZwl has brightened by a factor of ~2 (§ 2) in continuum flux and both the 
line strengths and the continuum slope have changed since the IUE data were taken, this is not a 
strong argument against the presence of Fe in in the HST spectrum. Due to the uncertainty in the 
identification and strength of this possible iron emission feature, it is not included in the Fe in and 
Fell templates. 

A1807-1875 emission: Fein (UV97,117 at AA1830-1855A and UV53,63 at AA1850-1871A; see 
Fig. 4) is the most likely identification, but fainter Fen (UV65,66) emission is also consistent 
with some of the residuals in the range AA1807-1875A (Fig. 5d). Fe II is, however, expected 
to dominate in the range AA1807-1836A (most of the weaker, blended iron multiplets are not 
marked in Fig. 4). Due to the heavy blending, no attempt were made to separate the individual 
Fen and Fein contributions. A compromise was made: the emission at AA 1837-1872A is treated 
entirely as Fe in emission (and included in the Fe in emission model and Fe in template) while the 
AA1802-1836A wavelength range is adopted as pure Fell. This is clearly an approximation. 

Al in AA1854, 1863 emission: The doublet is relatively strong, and well resolved in this spectrum 
(Fig. 5d). A broad base component is clearly present along with a number of iron emission transi- 
tions. We detect no obvious excess blueshifted emission and find that the most likely identification 
of the non-Aim emission in the A1850— 1875A region is that of iron (Figs. 4 and 5d). 

Si in] A1883 emission: This was identified by L97 as possible Sim] emission (Fig. 5d). The 
feature is observed at A1880.3A. The Fe n UV126 multiplet is also expected (and consistent with the 
emission) at AA1864. 6, 1864. 7,1880. 97, providing an alternative identification. The 1864A lines are 
blended severely with the Al ill emission and the identification cannot be confirmed via line strength 
arguments. When in doubt we prefer to underestimate the iron emission (see the introduction to 
this section and § 3.3.2), so we do not include this feature in the template. 

Fein A1892-1906 emission: The residuals between A1898 and 1907 were not fitted with Gaussian 
components (Fig. 5d). However, they were isolated by subtraction of fitted non-iron lines and then 
included in the Fein template (Fig. 3). With the Fein UV34 A 1914 emission feature observed 
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at ~1912A, the 1895A transition is expected at 1893A, the position of a weak feature (Fig. 4). 
Given the apparent weakness of the feature at A 1893A the Fe in UV34 multiplet is not likely to 
dominate the strength of Sim] A1892, detected at ~1890A (Fig. 5e). The UV34 triplet (1895A, 
1914A, 1926A) may have optically thin line ratios 0.9:1.0:0.3 (§ 5.2) while all transitions have 
equal strengths in the optically thick limit. As discussed in section 5.2, our modeling of the Fein 
UV34 triplet show that no tight constraints can be placed on the relative line ratios. However, 
this modeling does indicate an optically thin ratio (cf. models B and C in Fig. 5e). The spectrum 
indicates that this triplet is emitted in the optically thin region given the apparent relative weakness 
of the 1895A and 1926A features. The relatively strong Siiii]/Ciii] line ratio combined with the 
relatively strong Fein emission indicates rather high BLR densities. 

C in] A1909 emission complex: The emission feature is a complex blend of emission lines, but 
the spectral resolution and the narrow widths of the line cores permit a separation of the narrow- line 
cores of the Sim] A1892, Cm] A1909 and Fein AA1914, 1926 lines (Fig. 5d). Two broad components 
of Sim] and Cm] could be fitted (Table 3), although their individual fits, including their strengths, 
are not well constrained. A number of other Fe m transitions are detected in and around this line 
complex (Fig. 4 and Table 4). 

FemA1914 emission: The emission at AA1907— 1918 is consistent with contributions from a 
number of intermediate to strong Fem multiplets; among them, the strongest are UV34, 101, and 
83 and the weaker ones are UV57, 135, and 108. This heavy blending is probably responsible for 
the residual UV34 (A1914) feature appearing stronger than, e. g., UV68 A1952, and UV50 AA1987 
- 1996 (Fig. 4). If UV51 (AA1915.1, 1922.8,1930.4,1937.3,1943.4) is present, it must be faint. See 
section 5.2 for discussion of the relative Fe m UV34 contribution. 

Nil] A2141 emission: The pseudo-continuum was used as the local continuum level (Fig. 2). 

A2160-2190 region: The absolute strengths of the emission and absorption lines are not straight- 
forward to determine in the ~2100— 2250A region due to the presence of dust absorption features 
(§ 2) and the uncertainty in the placement of the absolute continuum level (Fig. 2). According to 
the line lists of Nahar (1995), the expected Fen transitions are faint between 2100 and 2200A. 

AA2250-2290 region: The possible Fem UV73,153 (AA2258.1, 2274.7, 2277.6, 2285.7) emission is 
not fitted and not included in the Fe m emission model as the identification is uncertain, due to 
poor wavelength coincidence. 

Fe in UV47 AA 2418,2438 emission: The profile was fitted with three Gaussian components, 
one of which accounts for a broad component and the weaker of the two narrow components is 
blueshifted ~500kms~ 1 with respect to the stronger one (Fig. 5f). The choice of including a broad 
component was based on the appearance of the Ha profile (Phillips 1977; L97). 
Unfortunately, the spectrum does not offer good constaints on the absolute strengths of the broad 
Fein UV47 component fit. A good fit requires that the actual continuum level for the fitting be 
placed ~6% below the observed flux level at the apparent line base. 
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The profile is also consistent with a fit of two narrow Gaussian components (and no broad com- 
ponent) and the actual continuum level (for the fitting) is placed at the observed flux level, but 
with the known similarity between the (Balmer) hydrogen lines and the Fen lines, the presence of 
a broad Fe in component is much more likely (see also L97) . 

A2435 emission feature: The identification is uncertain. Possibilities include the (slightly 
blueshifted) On (UV18) AA2433.6, 2444.4, 2445.6. Its blueshift (^Skins^ 1 ) is consistent with 
that of [Oll]A2470. Some Sii transitions are also expected nearby (Sii UV45 and UV2; Table 3). 
The feature was fitted with the pseudo-continuum as the local continuum. The best fit was ob- 
tained with two Gaussians (Fig. 5f); it is not clear how the individual Gaussian parameters (width 
and position) relate to those of the multiplet components. 

[On] A2470 emission: None of the lists of Fe n and Fe m transitions (including that by Giridhar 
& Arellano Ferro 1995) predict strong transitions at this wavelength. L97 identify the feature as 
[On] (Figs. 2 and 5f). The pseudo-continuum acts as the local continuum level in the component 
fitting. 

A2481 emission feature: The identification is uncertain. No obvious strong iron emission is 
expected at this position. C I A 2478 is a possible identification. The fit is based on the pseudo- 
continuum (Fig. 5f). 

Mg II AA2795, 2803 emission: The systemic redshift was measured with a two-Gaussian compo- 
nent fit (one doublet; § 3.3.2) to the narrow line core early in the course of this work. A more 
detailed and improved line fit reveals a second doublet with slightly larger line widths blueshifted 
400 km s -1 relative to the stronger doublet at rest in the quasar frame (Fig. 5g and Table 3). Each 
doublet has a thermalized line ratio (=1:1). 

The fit displays evidence for excess blueshifted emission similar to that in Ly a , N v, He n and 
possibly Siiv+Oiv], and Civ. L97 fitted the entire doublet (broad and narrow components) using 
Ha as a template profile. Two such profiles with relative strengths 1.2:1 fitted the line well, and 
they argued for higher electron densities in the Mgn line gas based on the thermalized doublet 
ratio (~1:1). 

It is clear from Figure 5g that no iron emission is left in the A2770-2820A range after subtrac- 
tion of the Mgn fit. Fen emission is not expected to be strong at the position of Mgn (D. Verner, 
1997, private communication), but is not entirely absent either (see Fig. 4 and Verner et al. 1999). 
As the data cannot constrain the strength of the broad component well, we prefer to underestimate 
the Fe n strength to prevent an overcorrection of the iron emission when the iron template is applied 
(see § 3.3.2). This choice of the local underlying continuum is consistent with that found by L97. 



B. The 3C273 HST Data 



In § 4.3 HST archival data, covering the full UV range from Ly a to beyond Mgll, of the nearby 
radio-loud quasar, 3C273, were used to demonstrate the applicability of the IZwl based Fen and 
Fein templates to the iron emission in other AGNs. The 3C273 UV spectrum was generated by 
combining data from the FOS and the Space Telescope Imaging Spectrograph (STIS). We used 
the 1" aperture G130H, G190H, and G270H FOS grating spectra from 1991 January 16, 14, and 
15, respectively, which only partially cover the Fen bump around Mgn. The spectral range was 
extended longward of Mgn by combining these spectra with a 52x2 aperture G430L STIS grating 
spectrum observed on 1999 January 31. The latter spectrum was multiplied by a factor 1.142 to 
coincide with the average flux level of the FOS spectra. The G270H spectrum was truncated at 
3200A (observed) before being combined with the G430L spectrum, which then provided the Mg n 
line profile. There is a ~50A gap between the G130H and G190H spectra. An interpolation was 
performed using the average continuum flux levels in the ~40A ends of each spectrum, avoiding 
absorption lines, to facilitate the fitting of the iron emission. The spectra were resampled to a 
common dispersion of 2.37A/pix in the rest frame based on the G430L spectrum. The lower 
resolution of the G430L spectrum did not affect the iron fitting in any way. This was confirmed by a 
simultaneous fitting to the higher resolution spectrum made from the FOS spectra only (0.44A/pix, 
rest frame; no fitting to the full 'small iron bump' was possible, however) which showed no significant 
differences in the fitted models. This is most likely due to the relative broadness (~4000kms~ 1 ) of 
the lines in the spectrum (§ 4.3). 
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Fig. 1. — The non-iron emission model (solid line) overplotted on the continuum-subtracted IZwl 
spectrum (dotted line). The assumed continuum of the HST spectrum is a broken power-law, 
Fj, ~ v~ a with a break at A1716A (see text for details). The absorption lines are marked using 
letters a-z below the spectrum and listed in Table 2. The detected emission lines are labeled. 
Geocoronal Lya is seen at A1145A (restframe). 

Fig. 2. — The iron template (solid line; sum of Fe n and Fe in) overplotted on the I Zw 1 spectrum 
(dotted line). The detected emission lines are labeled. 

Fig. 3. — The Fein template spectrum (dotted line) overplotted on the Fen template spectrum 
(solid line). Most of the Fein features are fitted with Gaussian profiles. 

Fig. 4. — Suggested identifications of individual, mostly unblended UV Fen and Fein emission 
multiplets based on the positions and (rough) strengths listed by Moore (1950). Each label consists 
of the ionizatation level and the multiplet number separated by a hyphen. Multiplets with labels 
in parentheses are slightly blended. Square brackets denote multiplets which are suggested present 
by the fitting process (see Appendix A; '2-10' is expected only). Heavily blended multiplets are not 
marked. Multiplets with labels ending with an 'm' have relative transition strengths inconsistent 
with the optically thin laboratory strengths. Note that the positions in the spectrum of the multiplet 
transitions are sometimes blueshifted by 1-2 A relative to the laboratory wavelengths marked. This 
is particularily evident for FemUV34 at ~1912A (labeled '3-34m'). See § 5.2 and Appendix A 
regarding the relative strengths in this multiplet. 

Fig. 5. — Examples of line fitting. Each individual component is shown as well as their sum (dot 
- dashed line) superposed on the data (solid line). The non-iron emission is dotted while the 
fitted iron emission is dashed. Note, the fit is sometimes so good it entirely coincides with the 
spectrum, a) Ly a complex; note the strongly blue-asymmetric N v profile. This is also seen in the 
non-deblended Civ and Mgn profiles (App. A), b) Siiv+Oiv] A 1400 blend, c) Civ complex; '?': 
the identification of this feature is uncertain; see text, d) Fen UV191, and Cm] complex [Note, 
the Znii ISM absorption at 1910A. Both the sum of the non-iron emission line fits and the sum of 
non-iron and iron fits are shown in dot-dashed curves. The latter (full) sum coincides so well with 
the data in many places that it is not easily seen, e. g., around Sim], Cm], and Fein redward of 
Cm]. The fit to Fen UV 191 also coincides well with the data]. The fit to Fein UV34 is model 
B. e) Other sample model fits to Fein UV34 as discussed in text. The individual fits to Sim] and 
Cm] are not shown, but the residuals are. f) Fein UV47 A 2418, 2438, and the non-iron features 
nearby, and g) Mgn. 

Fig. 6.— The HST spectrum of 3C273 (Q1226+023) and the best fitted iron 'models'. The upper, 
large panel in (a) shows the original spectrum (shifted for clarity; see § 4.3) with the continuum 
overplotted (top) compared to the residuals (middle) after subtracting the best fitted iron 'model' 
(bottom). The latter is displayed separately below the large panel for clarity. Panel (b) shows the 
original spectrum, two iron 'model' fits, and the corresponding spectrum residuals. Neither of the 
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iron 'models' fits the ~2500A multiplets very well, but they fit reasonably well to the blue and red 
part of the "multiplet" , respectively Panel (c) shows these iron 'models' more clearly. See § 4.3 
for more details. The flux density is in units of 10 -13 ergs cm -2 s _1 A -1 . 

Fig. 7. — Sample fitting of the iron emission in typical high-redshift quasar spectra. In the upper, 
large panels the original spectrum (top) with the power-law continuum fit, is compared to the 
residual spectrum (continuum superposed; middle) after subtracting the best fit iron emission 
model (bottom). In the lower, smaller panels the iron models are plotted separately for visibility. 
Q0020+022 and Q0252+016 display pronounced 'pseudo-continuum' in the AA1500 - 1900 A range. 
The data are from M. Vestergaard et al. (2001, in preparation) with flux density in units of 10~ 16 
ergscm~ 2 s _1 A -1 . See § 4.3 for continuum parameters and details on the relative spectrum shifts 
valid for each quasar. 
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Fig. 1. — The non-iron emission model (solid line) overplotted on the continuum-subtracted IZwl 
spectrum (dotted line). The assumed continuum of the HST spectrum is a broken power-law, 
Fj, ~ v~ a with a break at A1716A (see text for details). The absorption lines are marked using 
letters a-z below the spectrum and listed in Table 2. The detected emission lines are labeled. 
Geocoronal Lya is seen at A1145A (restframe). 
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Fig. 2. — The iron template (solid line; sum of Fen and Fein) overplotted on the IZwl spectrum 
(dotted line). The detected emission lines are labeled. 
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Fig. 3. — The Fein template spectrum (dotted line) overplotted on the Fell template spectrum 
(solid line). Most of the Fein features are fitted with Gaussian profiles. 
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Fig. 4. — Suggested identifications of individual, mostly unblended UV Fell and Fein emission 
multiplets based on the positions and (rough) strengths listed by Moore (1950). Each label consists 
of the ionizatation level and the multiplet number separated by a hyphen. Multiplets with labels 
in parentheses are slightly blended. Square brackets denote multiplets which are suggested present 
by the fitting process (see Appendix A; '2-10' is expected only). Heavily blended multiplets are not 
marked. Multiplets with labels ending with an 'm' have relative transition strengths inconsistent 
with the optically thin laboratory strengths. Note that the positions in the spectrum of the multiplet 
transitions are sometimes blueshifted by 1-2 A relative to the laboratory wavelengths marked. This 
is particularily evident for FeinUV34 at ~1912A (labeled '3-34m'). See § 5.2 and Appendix A 
regarding the relative strengths in this multiplet. 
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Fig. 4. — continued 
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Fig. 5. — Examples of line fitting. Each individual component is shown as well as their sum (dot 
- dashed line) superposed on the data (solid line). The non-iron emission is dotted while the 
fitted iron emission is dashed. Note, the fit is sometimes so good it entirely coincides with the 
spectrum, a) Ly a complex; note the strongly blue-asymmetric N v profile. This is also seen in the 
non-deblended Civ and Mgn profiles (App. A), b) Siiv+Oiv] A 1400 blend, c) Civ complex; '?': 
the identification of this feature is uncertain; see text, d) Fen UV191, and Cm] complex [Note, 
the Znii ISM absorption at 1910A. Both the sum of the non-iron emission line fits and the sum of 
non-iron and iron fits are shown in dot-dashed curves. The latter (full) sum coincides so well with 
the data in many places that it is not easily seen, e. g., around Sim], Cm], and Fein redward of 
Cm]. The fit to Fen UV 191 also coincides well with the data]. The fit to Fein UV34 is model 
B. e) Other sample model fits to Fein UV34 as discussed in text. The individual fits to Sim] and 
Cm] are not shown, but the residuals are. f) Fein UV47 A 2418, 2438, and the non-iron features 
nearby, and g) Mgn. 
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Fig. 5. — (See caption above) 
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Fig. 6.— The HST spectrum of 3C273 (Q1226+023) and the best fitted iron 'models'. The upper, 
large panel in (a) shows the original spectrum (shifted for clarity; see § 4.3) with the continuum 
overplotted (top) compared to the residuals (middle) after subtracting the best fitted iron 'model' 
(bottom). The latter is displayed separately below the large panel for clarity. Panel (b) shows the 
original spectrum, two iron 'model' fits, and the corresponding spectrum residuals. Neither of the 
iron 'models' fits the ~2500A multiplets very well, but they fit reasonably well to the blue and red 
part of the "multiplet", respectively. Panel (c) shows these iron 'models' more clearly. See § 4.3 
for more details. The flux density is in units of 10 -13 ergs cm -2 s _1 A -1 . 
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Fig. 7. — Sample fitting of the iron emission in typical high-redshift quasar spectra. In the upper, 
large panels the original spectrum (top) with the power-law continuum fit, is compared to the 
residual spectrum (continuum superposed; middle) after subtracting the best fit iron emission 
model (bottom). In the lower, smaller panels the iron models are plotted separately for visibility. 
Q0020+022 and Q0252+016 display pronounced 'pseudo-continuum' in the AA1500 - 1900 A range. 
The data are from M. Vestergaard et al. (2001, in preparation) with flux density in units of 10~ 16 
ergs cm -2 s^ 1 A^ 1 . See § 4.3 for continuum parameters and details on the relative spectrum shifts 
valid for each quasar. 
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Fig. 7. — (See caption above) 
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Fig. 7. — continued 
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Table 1. Log of Observations of I Zw 1 (PG 0050+124) 



Grating 


Date 


A-range 


AA 


Res a 


Detector 


Aperture 


Exp. Time 


A Offsets 15 




(UT) 


(A) 


(A/pix) 


(A) 




(arcsec) 


(sec) 


(A) 


G130H 


02/13/94 


1087-1606 


0.251 


0.96 


blue 


0.86 


29 700 


0.36230 


G190H 


09/14/94 


1572-2312 


0.359 


1.39 


amber 


0.86 


6 030 


1.10670 


G270H 


09/14/94 


2222-3277 


0.511 


1.97 


amber 


0.86 


2100 


1.03755 



a Spectral Resolution 

b Wavelength offset determined from ISM absorption lines; as described in § 2 
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Table 2. Absorption Line Identifications 3 



Apparent 13 


Ion 


Aio6 C 


\ meas 
A obs 


\meas 
^rest 


Av d 


Flux e 


EW 


FWHM 


Mult 


Comments 


Feature 




(A) 


(A) 


(A) 


(km/s) 




(A) 


(km/s) 






& Position 






















A1109 (a) 


Cm* 


1175. 5[1] 


1177.6 


1109.8 


+523 


1.57 


0.4 


125 


UV4 




A1124 (bj 


ai ii 


-I -I r\r\ A rol 

1190. 4[2J 


1190.3 


1121.8 


—25 


2.17 


0.6 


260 


UV5 






Si ii 


"i "i no o rol 

1193. 3[2J 


1192.7 


1124.1 


— 150 


1.50 


0.4 


260 


UV5 






Si ii 


1194. 5[1] 


1193.6 


1125.0 


—226 


1.09 


0.3 


260 


UV5 






SI II 


1197. 4[1] 


1195.3 


1126.6 


—526 


1.02 


0.3 


260 


UV5 




A1135 (c) 


Ni 


1200.0 f [l] 


1200.2 


1131.2 


+50 


2.90 


0.8 


595 


UV1 


A1199.5, 1200. 2, 1200.8 




Mnll 


1203.8[1] 


1203.9 


1134.7 


+25 


1.58 


0.4 


570 


UV26? 






Si g 


1205.6[1] 






-420 










Alt. ID 




Sim 


1206. 5 [2] 


1205.7 


1137.3 


-200 


3.59 


1.0 


580 


UV2 




A1145 (d) 


Mnll s 


1 ni /i /" r-i i 
121U.b[lJ 


1210.2 


1140.6 


— 105 


4.70 


1.3 


770 








Mnll g 


1213.3[1] 


1213.0 


1143.3 


-65 


2.51 


0.7 


725 








Si g 


1214.3[1] 


1214.0 


1144.2 


-75 


6.80 


1.9 


1025 








Lya s 


1215. 7[2] 


1215.5 


1145.6 


-50 


0.016 


0.005 


135 


UV1 






Lya 6 


1215. 7[2] 


1217 A 


1147.4 


+417 


7.20 


2.0 


990 


UV1 


z =0.0014 




T S 

Lya B 


1215. 7[2J 


1220.5 


1150.3 


+1177 


6.40 


1.8 


1065 


UV1 


z =0.004 




Lya 8 


1215. 7[2] 


1225.7 


1155.2 


+2460 


1.45 


0.4 


725 


UV1 


z =0.0082 


A1187 (e) 


C IT 

a ii 


i o r f\ v rol 

1259. 5[2J 


1260.1 


1187.7 


+155 


2.69 


0.8 


325 


UV9 




A1226 (g) 


Oi 


1302. 2 [2] 


1302.1 


1227.2 


-23 


4.39 


1.2 


290 


UV2 






Si ii 


i on a a Tol 

1304. 4[2J 


1304.5 


1229.5 


+23 


4.23 


1.2 


290 


UV3 






Nv h 


1238.8[2] 


1306.0 


1230.9 


-1913 


5.23 


1.5 


290 


UV1 


z aba = 0.0064 


A1234 (n) 


Nv h 


1242. 8[2J 


1310.2 


1234.9 


— 1907 


6.15 


1.7 


340 


UV1 




A1257 (i) 


Cn 


1334.5[2] 


1334.3 


1257.6 


-42 


0.23 


0.06 


210 


UV1 






ClI* 


1335. 7[2] 


1334.5 


1257.8 


-264 


2.11 


0.6 


210 


UV1 




A1317 (j) 


Si iv 


1393.8[2] 


1394.0 


1313.9 


+43 


0.88 


0.3 


310 


UV1 






Si iv 


1402.8[2] 


1403.2 


1322.5 


+85 


0.38 


0.1 


310 


UV1 




\ 1 A A(\ C\r\ 

A144U [K) 


51 II 


1 C^Ofi 7T01 
lOZD. / [zj 


lOzD.o 


1 A QO n 

14oy.U 


+iy 


z.zz 


U.D 




T J\TO 

U V z 






Sin 


1533.4[1] 


1533.6 


1445.8 


+39 


0.92 


0.3 


495 


UV2 




A1460 (l) 


Civ 


1548. 2 [2] 


1548.0 


1459.0 


-39 


2.35 


0.7 


570 


UV1 






Civ 


1550.8[2] 


1551.2 


1462.0 


+74 


1.59 


0.5 


570 


UV1 




A1575 (m) 


Aln 


1670.8[2] 


1670.8 


1574.7 





1.08 


0.3 


245 


UV2 




A1703 (n) 


Sin 


1808.0[2] 


1807.9 


1704.0 


-15 


0.45 


0.1 


220 


UV1 




A1823 (o) 


Fei 


1934.5[1] 


1934.1 


1822.9 


-60 


0.53 


0.2 


220 


UV37 






Fei 


1937.3[1] 


1936.4 


1825.1 


-140 


0.49 


0.2 


2190 


UV35 
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Table 2 — Continued 



Apparent 13 


Ion 


Ala6 C 


\ meas 
A obs 


\meas 
surest 


Av d 


Flux e 


EW 


FWHM 


Mult 


Comments 


Feature 




(A) 


(A) 


(A) 


(km/s) 




(A) 


(km/s) 






& Position 






















A1909 (p) 


Znn 


2026. 1[2] 


2026.1 


1909.6 





0.97 


0.3 


180 


UV1 




A2121 (q) 


Fen 


2249.9[1] 


2248.8 


2119.5 


-145 


0.41 


0.1 


305 


UV5 






Fel 


2251. 5[1] 


2250.7 


2121.3 


-105 


0.5 


0.2 


265 


UV16 






Nil g 


2255 n 1 


2253.7 


2124.1 


—240 


0.23 


0.1 


295 


UV34 






Ni I 


2256 fifll 


2256.9 


2127.1 




0.23 


0.1 


280 


UV9 






Fel 


9960 9M1 






—400 








UV15 


Alt. ID 




Fe II 


99fi0 R\9] 

.0 1 ^1 


2260.7 


2130.7 


— 15 


0.22 


0.1 


285 


UV4 




A2171 (r) 


Mnll g 


2299. 7[1] 


2300.9 


2168.6 


+155 


0.23 


0.1 


280 


UV2 






Mnn g 


2305. 7[1] 


2304.5 


2172.0 


-155 


0.42 


0.2 


280 


UV2 




A2206 (s) 


Nil 


2338. 2[1] 


2339.5 


2205.0 


+165 


0.51 


0.2 


290 


UV8 






Fen 


2344.2[2] 


2344.0 


2209.2 


-30 


1.34 


0.5 


290 


UV3 




A2235 (t) 


Coi g 


2365.8[1] 


2365.4 


2229.4 


-50 


0.29 


0.1 


295 


UV6 






Crl 


2365 5M1 






— 10 










Alt. ID 




All 


2367.8[1] 


2368.3 


2232.1 


+60 


0.48 


0.2 


295 


UV4 






All 


2373.9 f [l] 


2372.0 


2235.6 


-240 


0.26 


0.1 


295 


UV4 






Ah 


2372 8M1 






— 100 








UV3 


Alt. ID 




Fen 


2374.5[2] 


2374.5 


2238.0 





0.92 


0.4 


295 


UV2 




A2245 (u) 


Fen 


2382.8[2] 


2382.9 


2245.9 


+ 10 


1.49 


0.6 


230 


UV2 






Coi 


2385.6[1] 


2386.0 


2248.8 


+50 


0.2 


0.1 


145 ! 


UV6 




A2268 (v) 


Coi g 


2408.0[1] 


2407.6 


2269.2 


-50 


0.48 


0.2 


245 ! 


UV6 




A2445 (w) 


Fen 


2586. 7[2] 


2586.9 


2438.2 


-25 


0.4 


0.2 


155 ! 


UV1 






Fen 


2600. 2 [2] 


2600.1 


2450.6 


-10 


1.84 


0.8 


285' 


UV1 




A2635 (x,y) 


Mgn 


2796.4[2] 


2796.4 


2635.6 





2.15 


1.0 


215 


UV1 






Mgn 


2803. 5 [2] 


2803.7 


2642.5 


+20 


2.14 


1.0 


215 


UV1 




A2688 (z) 


Mgi 


2853.0[2] 


2853.1 


2689.1 


+ 15 


0.61 


0.3 


215 


UV1 





a All the absorption is consistent with Galactic (ISM) absorption, except where noted; see "h". 
Cm*and Cn* are meta-stable transitions. 

b A simultaneous fit was made to each line (in column 2) listed as a Feature (in column 1) unless otherwise noted. 
Note, the wavelength of each Feature refers to the wavelength as it appears in the restframe spectrum in Figure 1. 

References. — [1] Morton 1991; [2] Savage et al. 1993. 

d Av is the line peak velocity shift from the expected rest frame position. The resolution is 470 km s _1 at A1200 A, 355 km s _1 
at A1600A, 280 kms -1 at A2000 A, 240 kms -1 at A2400 A and 200kms _1 at A2800 A (Leitherer 1995). 

e Units of 10~ 14 ergs cm -2 s _1 . 

f Multiplet wavelength or mean wavelength of more than one transition. If more than two transitions, the wavelengths will 
be listed in the Comments column. 
s The line identification is uncertain. 

h Not Galactic absorption. This is absorption intrinsic to the object. Laor et al. 1997b identified these absorption features 
as UV (Nv) Quasar associated absorption. 

'A simultaneous fit was not possible due to the relative locations of the absorption lines. 
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Table 3. Emission Line Fits and Measurements 



Ion 


Mab 


\ meas 
surest 


A b 


Line b lux 


EW d 


b WrlM 


Comments 




(A) 


(A) 


(kms^ 1 ) 




(A) 


(kms^ 1 ) 




Cm* 


1175.5 


1173.9 


-408 


9.46 


2.7 


1265 




Sim 


1206.5 


1205.9 


-150 


11.89 


3.4 


525 




' * total 


1215.7 






694.99 


197.4 






— narrow 


1215.7 


1210.8 


-1209 


69.19 


19.6 


1020 






1215.7 


1214.4 


-320 


182.00 


51.7 


1025 




— broad 


1215.7 


1216.9 


+296 


443.80 


126.1 


11025 




Ov ? 


1218.3 


1219.0 


+172 


2.94 


0.8 


395 


ref[l], ID uncertain, see App. A 




1240.2 






89.66 


25.50 






Nv 


1240. T 


1230.1 


-2443 


17.84 


5.1 


1200 


A1238.8,1242.8 




1238.8 


1234.9 


-944 


37.79 


10.7 


1200 






1242.8 


1239.5 


-797 


34.03 


9.7 


1200 




Si Htotal 


1249 






4.48 


1.2 








1248.4 


1249.7 


+310 


1.50 


0.4 


835 


Sgg AppGiidix A 




1251.1 


1250.7 


-95 


2.98 


0.8 


835 




Si Utotal 


1263.3 






26.92 


7.7 








1260.4 


1257.2 


-760 


4.16 


1.2 


775 






1264.7 


1260.7 


-950 


10.47 


3.0 


860 






1265.0 


1264.2 


— 190 


12.29 


3.5 


845 




Si III tota ; 


1297 






4.54 


1.3 




At)=-2000kms _1 if Sill+Ol 


Sim* 


1295.5 


1293.9 


-370 


2.56 


0.7 


595 


ID uncertain, See App. A 




1298.9 


1296.5 


—555 


1.98 


0.6 


590 




Si 11+ Oltotal 


1306 






27.92 


8.1 






Si n+ Oi 


1302.2 


1299.0 


-735 


3.04 


0.9 


755 






1304. 6 e 


1302.1 


-575 


8.67 


2.5 


725 


A1304.4,1304.9 




1306.0 


1305.5 


-115 


10.37 


3.0 


710 






1309.3 


1309.3 





5.84 


1.7 


820 




C Utotal 


1335.3 






12.53 


3.6 






Cn 


1334.5 


1330.8 


-830 


3.39 


1.0 


890 






1335.7 


1334.0 


-380 


9.14 


2.6 


890 




A1343 ID? 




1343.9 




3.14 


0.9 


955 


Possible IDs below: 
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Table 3 — Continued 



ion 


*lab 


\meas 
^rest 


Z_i u 


T itwi T7 1 1n-v < -' 
Ijlllfc! r 1LIX 


r. vv 


r vv jnivi 


L> omment s 




\ A ) 








\ A ) 


^Klll S J 




Oiv 


1343.5 


1343.9 


-90 








ref[l] 


Can 


1342. 3 e 


1343.9 


+360 








A1341. 9,1342. 6 


Fein 


1343. 2 e 


1343.9 


+ 155 








A1338,42. 5,43. 3,45. 7,46. 4,48.0; r 


Fev 


1343.1 


1343.9 


+ 180 








Alt. ID 


Fev 


1345.6 


1343.9 


-380 








Alt. ID 


Oi 


1355.6 


1343.9 


-2590 








Alt. ID 


SiW+Ow]total 


1400.0 






96.74 


27.8 






— narrow 


1395. 5 e 


1387.4 


-1740 


5.67 


1.6 


920 


A1393. 8,1397.2 


— narrow 


1399.8 


1392.2 


-1630 


10.57 


3.0 


1020 




— narrow 


1401. 8 e 


1397.5 


-920 


11.99 


3.5 


1010 


A1401. 2,1402.8,1404.8 




1407.4 


1401.2 


-1320 


16.73 


4.8 


1010 




— broad 


1400. 5 e 


1393.4 


-1520 


51.78 


14.9 


7760 


A1399. 8,1401. 2 


A1520 


1526.7 


1521.2 


-1080 


1.73 


0.5 


478 


ID? 




1533.4 


1524.0 


-1840 


0.86 


0.3 


478 


Shifts as if Si n 


Si Utotal 


1531.2 






23.21 


6.7 






Sin 


1526. 7 f 


1526.3 


-75 


8.26 


2.4 


1120 






1533.4 f 


1532.0 


-270 


14.95 


4.3 


1120 




A154CW 


1549.1 






11.56 


3.3 




ID? 


A1540 


1548.2 


1535.8 


-2400 


6.65 


1.9 


630 


Shift as if C iv 




1550.8 


1539.3 


-2225 


4.91 


1.4 


630 




C IVtotal 


1549.1 






74.68 


21.9 






C IV narrow 


1548. 2 f 


1543.3 


-950 


29.11 


8.6 


1100 


Best fit 


— narrow 


1550. 8 f 


1548.4 


-465 


14.93 


4.4 


895 




— broad 


1549.1 


1538.1 


-2130 


30.64 


8.9 


12265 




C IV narrow 


1548. 2 f 


1543.3 


-950 


27.73 


8.1 


1070 


See App. A 


— (alternative fits) 


1550. 8 f 


1548.4 


-465 


14.10 


4.1 


1070 








1547.9 




2.29 


0.7 


535 


Fit to residual 


C IV complex 8 


~1548 






126.78 


36.47 




A1505-1560A 
















(incl. iron) emission 


Hell tota i 


1640.4 






11.24 


3.3 






Hen 


1640.4 a 


1628.4 


-2195 


5.38 


1.6 


840 






1640.4 a 


1633.1 


-1335 


3.92 


1.1 


840 
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Table 3 — Continued 



Ion a 




\meas 
surest 




Line Flux c 


EW d 


FWHM 


Comments 




(A) 


(A) 


(krns" 1 ) 




(A) 


(krns" 1 ) 




Hen 


1640. 4 f 


1639.2 


—220 


1.94 


0.6 


820 




Om] 


1664. T 


1663.7 


-90 


2.61 


0.8 


940 


Fit with Al II 


Alii 


1670. 8 a 


1669.7 


-200 


2.15 


0.5 


940 


Fit with O in] 


Niu] tota i 


1750.5 






16.85 


4.9 






Niii] 


1746.8 


1737.0 


-1685 


2.47 


0.7 


830 






1748.7 


1741.0 


-1320 


3.16 


0.9 


830 






1749.7 


1744.0 


-975 


2.82 


0.8 


830 






1752.2 


1746.8 


-925 


5.25 


1.6 


830 






1754.0 


1751.0 


—575 


3.15 


0.9 


830 




& l ll total 


1814.0 






4.98 


1.5 






Sin 


1808.0 f 


1807.6 


-65 


2.06 


0.6 


885 






1817.2 f 


1816.1 


-180 


2.92 


0.9 


925 




lll total 


1857 






30.10 


9.6 






— narrow 


1854.4 


1852.7 


-275 


3.50 


1.1 


870 




— narrow 


1862.8 


1861.3 


-240 


3.65 


1.2 


910 




— broad 


1854.4 


1849.1 


-860 


13.09 


4.2 


4530 






1862.8 


1856.0 


— 1095 


9.86 


3.1 


4510 




[Sim] 


1882.7 


1880.3 


-380 


1.02 


0.3 


340 


FenA1881?, Model B 


SilIl]total 


1892.0 






29.25 


9.6 




Model B 


11 al 1 \J W 


1892.0 


1889.5 


—400 


11.80 


3.9 


800 




— broad a 


1892.0 


1889.5 


-400 


17.45 


5.7 


3950 


Fit with C III] 


C lll]total 


1908.7 






51.85 


17.0 




Model B 


— narrow 


1908.7 


1906.8 


-300 


9.24 


3.0 


880 




— broad a 


1908.7 


1906.8 


-300 


42.61 


14.0 


5500 


Fit with Si in] 


C III] complex I 


-1909 






82.12 


26.9 




Ciii]+Siiii]+[Siin] 


C III] complex II 


-1909 






112.22 


36.5 




Ciii]+Siiii]+[Siin]+Aliii 


Nil] 


2141.4 


2141.0 


-55 


2.66 


1.0 


935 




C ll]total 


2326.6 






3.11 


1.2 






Cu] 


2325.4 


2324.0 


-180 


0.95 


0.4 


425 






2326.1 


2325.6 


-65 


1.05 


0.4 


425 
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Table 3 — Continued 



Ion a 




\meas 
^rest 


Av h 


Line Flux c 


EW d 


FWHM 


Comments 






{ A ) 


(Km s ) 






(km s ) 






2327.7 


2327.2 


-65 


0.55 


0.2 


425 






2328.8 


2328.7 


-10 


0.56 


0.2 


425 




A2435 ID? 




2435.4 




3.18 


1.34 




model sum 






2436.8 




1.40 


0.59 


430 


modeled 






2433.5 




1.78 


0.75 


430 


modeled 


On(UV18) 


2440.7 6 


2435.4 


-650 








Suggested ID,ref[3] 


Sii(UV45) 


2435.2 


2435.4 


+25 








Alt. ID,ref[3] 


Sii(UV2) 


2444.3 e 


2435.4 


-1090 








Alt. ID,ref[3] 


[On] 


2470.3 


2464.6 


-690 


3.05 


1.0 


665 


ref[l] 


Ci ? 


2478.6 


2482.0 


+410 


2.52 


1.1 


690 


ID uncertain,ref[3] 


A 1 TT 

Al II 


ZD / U.U 


ZOO / .0 


— Z40 


Q 1 A 
0.14 


1.0 


SOU 




Mgll totQl 


2797.9 e 






94.31 


45.8 




A2795.5,2802.7 


— narrow 


2795.5 


2796.2 


+75 h 


12.13 


5.9 


640 


4 Gauss fit 




2802.7 


2803.4 


+75 h 


12.13 


5.9 


640 


ditto; doublet at rest 




2795.5 


2792.0 


-375 


9.06 


4.4 


790 


ditto; shifted doublet 




2802.7 


2799.4 


-350 


9.06 


4.4 


790 


ditto 


— broad 


2797.9 e 


2793.1 


-515 


51.94 


25.2 


5720 


A2795.5,2802.7 


Mgll narrow 


2795.5 


2796.3 


+85 h 


33.84 


16.5 


1130 


2 Gauss fit; for z 


(alternative fit) 


2802.7 


2803.5 


+85 h 


8.98 


4.4 


550 


measurement 



a Each line/blend component, identified by element name in column 1, was fitted separately from the other (named) 
components - unless otherwise noted. Unnamed components following a named one were fitted simultaneously with 
this component. See also Appendix A for further details on each fitted feature. 

h Av is the line peak velocity shift from the expected rest frame position. The resolution is 470kms _1 at A1200 A, 
355 km s" 1 at A1600 A, 280 km s" 1 at A2000 A, 240 km s" 1 at A2400 A and 200 km s" 1 at A2800A (Leitherer 1995). 

c Units of 10~ 14 ergscm^s" 1 . 

d Some line fits may have required a local continuum level different from the global continuum setting, but the EW 
measurements are always with respect to the global continuum. 

e Multiplet wavelength or mean wavelength of more than one transition. 

f A separate fit was made to each component; non-simultaneous from the other transitions. 

g The C iv complex contains all the emission above the continuum from A1505 to 1560A: A1505 - 1525A (incl. iron) 
emission, SinAA1526, 1533, A1540A feature, and the broad and narrow components of Civ. 

h Insignificant velocity shift. This shift is due to the profile fitting uncertainty and is not a real effect. 

References. — Wilkes 2000, and Morton 1991 unless specifically marked: [1] Penston et al. 1983, [2] Kurucz 1997, 
[3] Moore 1950. 
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Table 4. Fitted Iron Emission Features a 



Feature 13 


\meas 
^rest 

(ky 


A/ i 

(A) 


Ion 


Mult 


Av d ' e 
(kms" 1 ) 


Flux f 


EW 

(A) 


FWHM e 
(kms" 1 ) 


Comments 


A1495 


1493.2 


1493.6 g 


Fein 


UV85 


-80 












1494.1 


1495.2 


Fein 





-220 


3.40 


1.0 


710 






1498.6 


1498.8 


Fein 


— 


-40 


2.87 


0.8 


710 




A1505 


1502.2 


1503.0 


Fein 




-160 


2.59 


0.8 


678 






1505.4 


1505.2 


Fein 


UV85 


+40 


1.10 


0.3 


678 






1509.2 


1509.2 


Fein 








1.20 


0.3 


678 




A1515 


1513.9 


1514.5 


Fein 




— 120 


4.85 


1.4 


900 






1518.5 


1518.8 


Fein 




-60 


5.00 


1.5 


836 




A1520 


1521.2 


1524.5 


t-i h 

Fein 




-650 


1.73 


0.5 


478 


Iron or Sim? 




1524.0 


1525.0 


Fein h 




-195 


0.86 


0.3 


478 








1525.8 


Fein 





-355 








Alt. ID 


A1786tota! 


1785.4 


1786.0 


Fen 


UV191 


-100 






550 


narrow, 1 Gaussian 


\1786total 


1783.0 


1785.3 


Fen 


UV191 


-385 


1.53 


0.5 


445 


3 Gauss model ref[4] 




1785.6 


1786.7 


Fen 


UV191 


-185 


2.63 


0.8 


445 


ditto 




1788.6 


1788.0 


Fen 


UV191 


+100 


0.85 


0.3 


445 


ditto 


A1850 


1850.3 


1850.0 




UV53,63 


+50 








Possibly blended w/ 






1849.8' 


Fein 


UV97 


+80 


0.45 


0.1 


265 


A1849.4,1850.2 


A1855 


1856.2 


1854.9' 


Fein 


UV63 


+210 


0.41 


0.1 


295 


A1854. 8, 1855.0 




1857.0 


1856.7 


Fein 


UV63 


+50 


0.29 


0.1 


265 






1858.2 


1858.5 


Fein 


UV63 


-50 


0.24 


0.1 


140 




A1870 


1867.9 


1866.5' 


Fein 


UV52 


+225 


0.20 


0.1 


260 


A1866.3,1866.6 






1867.6' 


Fen 




+50 








A1867.2, 1867.9 Alt. : 




1869.7 


1869.8 


Fein 


UV52 


-15 


0.22 


0.1 


140 








1869.6 


Fen 




+15 








Alt.ID 




1871.0 


1871.2 


Fein 


UV52 


-30 


0.09 


0.03 


95 








1871.0 


Fen 













Alt. ID 


A1885 


1884.6 


1885.6' 


Fein 


UV96 


-160 


1.808 


0.6 


640 


A1885. 1,1886.0 




1886.0 


1886.8 


Fein 


UV52 


-125 


1.43 


0.5 


640 






1887.1 


1887.5 8 


Fein 


UV52 j ;53 j 


-65 











A1895 1892.7 1894.8 Fern UV34 -330 

1895.5 Fern UV34 -445 

1894.0' Fem h UV96,83 -205 

1896.6 1897.1' Fern UV83 -80 

1896.8 UV83 -30 

1901.0 1901. l g Fern UV95 -15 



3.76 1.2 Model B ref[3] 
ref[m] 

■■■ A1892. 9, 1895. 0(96), 1894.0(83) 

w/Sim 1892 + UV96,83? 

0.60 0.2 310 A1896.8,1897.4(Kurucz) 

(Moore) 



A1912 1912.0 1914.1 Fern UV34 -330 

1911.3 Fern UV135 +110 



10.04 3.3 Model B See also App. A ref[4] 
■ ■ ■ Possibly blended w/ UV34,57 
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Table 4 — Continued 



Feature 13 


\meas 
vest 


A; a 6 


Ion 


Mult 


Av d,e 


Flux f 


EW 


FWHM e 


Comments 




(Ay 


(A) 






(kms -1 ) 




(A) 


(kms" 1 ) 








1912.9 


Fein 


UV57 


-140 








weak, blended w/ UV34,135 






1915.1 


Fein 


UV51 


-485 








Alt. ID ref[5] 






1916.5 


! Fein 


UV95 










blended 






1917.5 


Fein 


h,fc TTV1D1 

U V 1U1 


— ouu 








A If TT") 




1917.6 


1917.4 8 Fein 


UV95 


+30 












1922.7 


1923.0 


s Fein 


UV95 


-45 










A1924 


1923.6 


1923.5 


Fein 


UV95,57 


+15 


1.49 


0.5 


675 


A1923.0(95), 1923.9(57) 




1924.5 


1924.5 


Fein 


UV79 





3.36 


1.1 


750 






1923.9 


1926.3 


Fein 


UV34 


-375 


4.24 


1.4 


Model B 






1925.5 


1926.0 


Fen 


UV123 


-75 








blended w/UV34,57? 






1926.0 


' Fein 


UV57 










blended 






1928. 3 g Fein 


UV95 










blended 




1932.8 


1932.8 


s Fein 


UV95 













A1935 


1935.9 


1935.3 


Fen 


UV96 


-90 


0.76 


0.3 


385 


UV96:1936.8fl917.3 blended 




1938.7 


1938.8 


Fein 


UV95 


-15 


0.60 


0.2 


385 






1941.0 


1941.1 


Fein 


TTV7Q 

U V ID 






n 9 


GOO 


\ 1 Q40 fi 1 Q41 fi 


A1950 


1948.0 


1948.4 


Fen 


UV123 


-60 


2.40 


0.8 


830 






1 Q4Q 7 


1949. 7 g Fein 




n 
u 












1951.6 


1951.3 


Fein 


UV68 


+45 


2.39 


0.8 


825 


[m]:A1951.0, Av=+90 




1953.1 


1953.0 


Fein 


UV68 


+15 


1.70 


0.6 


825 


A1952.6,1953.3 






1953.3 


Fein 


TTV89 










hlpnHprl w/TTVRS? 




1953.5 


1953.5 


Fein 


UV82 















1957.9 


1957.9 


Fein 


UV147 





2.04 


0.7 


825 






1958.8 


1958.6 


! Fein 


UV55 


+30 














1960. 3 g Fein 


UV82 










blended 




1963.6 


1964.3 


Fein 


UV82 


-100 


1.67 


0.6 


820 


A1964.0,1964.2,1964.8 


A1970 


1970.1 


1969.9 


Fein 


— 


+30 


0.25 


0.1 


265 




A1975 


1973.5 


1976.1 


Fein 


h UV54 


-395 


0.35 


0.1 


255 


ID? 




1975.4 


1976.1 


Fein 


UV54 


-105 


0.42 


0.1 


255 




A1990 


1977.5 


1977.3 


_ 




+30 


0.14 


0.1 


190 






1981.5 


1982.5 


Fein 


UV54,56 


-150 


2.17 


0.7 


735 


A1982.1(54), 1982.8(56) 




1985.5 


1985.0 


Fein 


UV56,81 


+80 


1.00 


0.3 


730 


A1985.1,1987.8(57), 1984.3(81 




1989.6 


1989.5 


Fein 


UV50 


+15 


2.41 


0.8 


730 


A1987.5,1990.0,1991.6 




1994.2 


1994.2 


Fein 


UV50 





1.85 


0.6 


730 


A1993.3,1994.1, 1995.3 




1994.2 


1992.1 


Fein 


UV81 


+300 


1.85 


0.6 


730 


A1992. 0,1992. 2 




1999.2 


1999.8 


Fein 


UV55,81 


-90 


1.48 


0.5 


725 


A1999.6(55,81),2000.2(81) 






1999.4 


Fen 


UV186 


-30 








blended w/UV55,81? 




2004.5 


2004.1 


Fen 


UV83 


+60 


1.79 


0.6 


725 








2003.5 


Fein 


UV55 


+150 








blended w/UV83? 




2007.2 


2007.8 


' Fein 


UV55 


-90 














2017.9 


! Fen 


UV186 










blended 






2027.8 


s Fen 


UV186 










blended 
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Feature 13 


\meas 
^rest 

(Ay 


Aiat 
(A\ 


Ion 


Mult 


Av d ' e 
(km s _1 ) 


Flux f 


EW 

(A) 


FWHM e 
(kms -1 ) 


Comments 


A2030 


2032.5 


2033.1 






-90 


0.22 


0.1 


255 


ID? 


A2035 


2035.3 


2035.1 


Fen 


UV186 


+30 


0.97 


0.3 


410 


(Kurucz) 






2034.5 






+115 








(Moore) 




2036.6 


2036.8 s 


Fein 


UV60 


-30 










A2040 


2038.5 


2038.9 


Fein 


UV60 


-60 


0.83 


0.3 


410 








2039.5 


Fein 


UV134 


-145 








Alt. ID 




2044.7 


2045.0 s 


Fein 


UV60 


-45 












2049.8 


2050.7 s 


Fein 


UV60 


-130 










A2060 


2055.6 


2056.8 


Fein 


UV105 m 


-175 


2.25 


0.8 


815 


[m]: A2055.9, Av=-45 




2057.1 


2056.9 s 


Fein 


UV78 


+30 












2059.7 


2059.7 s 


Fein 


UV78 















2060.2 


2058.6 


Fein 


UV100 


+160 


1.98 


0.7 


810 


blended w/UV78 A2059.7? 




2061.3 


2061.6 


Fein 


UV48 


-45 


1.18 


0.4 


805 


blended w/UV78 A2061.8? 




2061.8 


2061.8 s 


Fein 


UV78 















2067.2 


2068.2 


Fein 


UV48 


-145 


3.59 


1.3 


810 






2077.6 


2077.8 s 


Fein 


UV105 


-30 










A2080 


2078.3 


2079.0 


Fein 


UV48 


-100 


3.11 


1.1 


735 


ref[2] 


A2100 


2083.0 


2083.1 


Fein 




-15 


0.42 


0.2 


280 








2083.5 


Fe III h 


UV124 


-70 


0.42 


0.2 


280 


Alt. ID 




2085.9 


2085.0 


Fein 


UV67 


+130 


0.79 


0.3 


375 






2086.0 


2086.1 s 


Fein 


UV105 


-15 












2088.8 


2088.7' 


Fein 


UV67,77 


+15 


0.62 


0.2 


320 


A2087.8, 2088. 6, 2089. 8(77), 2088. 






2088.1' 






+100 








[m]:A2087.1, 2087.9,2089. 1(77) 




2091.7 


2092.0 


Fein 


UV77 


-45 


0.31 


0.1 


305 


[m]:A2091.3, Ad=+60 




2094.7 


2094.2 


Fein 


UV77 


+70 


0.57 


0.2 


255 


[m]:A2093.5, Au=+170 




2095.0 


2095.3 s 


Fein 


UV105 


-40 












2097.9 


2098.3 s 


Fein 


UV67,66 


-60 


0.76 


0.3 


400 


A2098.2,2098.4 








Fein 


UV67,66 










[m]: A2097.5(67),2097.7(66) 




2100.7 


2100.8' 


Fein 


UV129 


-15 


0.25 


0.1 


280 


A2100.0,2101.6 








Fein 


UV129 










[m]: A2092.9 S 2099.3, 2101.0 




2103.4 


2103.8 


Fein 


UV66 j 


-60 


0.47 


0.2 


380 






2105.8 


2105.0 


Fein 


UV146 


+115 


0.85 


0.3 


350 








2107.3 


Fein 


UV66 


-215 








blended w/UV146? 




2108.4 


2108.7 


Fein 


UV105 


-40 


0.32 


0.1 


340 








2107.6 


Fen 


UV250 


+115 








blended w/Fem UV105? 




2109.6 


2109.4 


Fein 


UV105 


+30 


0.36 


0.1 


375 






2113.2 


2114.3' 


Fein h 




-155 


0.90 


0.3 


465 


A2114.0,2114.6 


A2418 


2418.4 


2418.6 


Fein 


UV47 


-25 






715 


narrow, 1 Gauss 


\241Stotal 




2418.6 


Fein 


UV47 




12.27 


5.1 




3 Gauss fit,ref[l] 


A2418 


2415.0 


2418.6 


Fein 


UV47 


-445 


1.20 


0.5 


590 


narrow 




2418.6 


2418.6 


Fein 


UV47 





3.61 


1.5 


475 


narrow 
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Table 4 — Continued 



Feature 13 


\ meas 
surest 




Ion 


Mult Av A ' e 


Flux f 


EW 


FWHM e Comments 




(Ay 


(A) 




(kms^ 1 ) 




(A) 


(kms^ 1 ) 




2418.8 


2418.6 


Fe in 


UV47 +25 


7.46 


3.1 


1910 broad 


A2438 


2439.6 


2438.9 


Fe in 


UV47 +85 


1.17 


0.5 


300 [m]: A2438.2, Au=+170 



a The features listed here are those which are strong enough to be fitted and for which fitting was required in the 
deblending process of the non-iron lines. See also Appendix A for details on each fitted feature. 

b A simultaneous fit was made to each line blend listed as a Feature (column 1) unless otherwise noted. 

c If no wavelength is listed the entry is an alternative identification to the entry above it. 

d Aw is the line peak velocity shift from the expected rest frame position. 

e The resolution is 470 km s" 1 at A1200A, 355 km s" 1 at A1600A, 280 km s" 1 at A2000 A, 240 km s" 1 at A2400A 
and 200 km s" 1 at A2800 A (Leitherer 1995). 

f Units of 10~ 14 ergs cm -2 s -1 . 

g The transition is detected only and/or is likely present; feature is not modeled and measured. 
h The line identification is uncertain. 

'Multiplet wavelength or mean wavelength of more than one transition. The wavelengths are listed in the Comments 
column. 

J Other transitions of this multiplet are possibly present, but (heavily) blended. 

k This transition may also contribute, but it is likely much weaker than the transition strength indicate. 
'A separate fit was made to this component; non-simultaneous from the other transitions. 
"Transitions without multiplet number contribute as well. 

References. — Moore 1950 ([m]); Kurucz 1997 (air-wavelength for A> 2000A), unless otherwise noted. [1] Ekberg 
1993, [2] Hartig & Baldwin 1986, [3] Marziani et al. 1996, [4] Penston et al. 1983, [5] Fuhr ct al. 1988. 
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Table 5. Pure Iron Emission Windows 



A limits a 


1 


2 


3 


4 


Range 
5 


number 
6 


*7 
I 


8 


9 


10 


11 


Ai (A) 
A 2 (A) 


1350 
1365 


1427 
1480 


1490 
1505 


1705 
1730 


1760 b 
1800 


1930 c or 1942 
2115 


2250 
2320 d 


2333 d 
2445° 


2470 
2625 


2675 
2755 


2855 
3010 



a The precise limits vary from object to object 

b Avoid Nm] A 1750 line. The width and position may vary from object to object 

c Avoid C in] A 1909 broad base if present. If this base is strong, use the A 1942 — 2115 A range. 

d Avoid C n] A2326 

e Exclude [Ne iv] A2419 if present 



- 73 - 



-74- 



Table 6. Model Fit Measurements of Fein UV34, [Sim], Sim], and Cm]. 



Model 


Line 


-Mab 


L/ine r lux/ 1U 


HjW 


r Wrilvl 


Flux Ratio 






( A ) 


(ergs cm - ^ s — * A - 




(km s — 




Model A 
















Fe m UV34 


1895.5 


4.55 


1.5 




1.0 






1914.1 


4.50 


1.5 




1.0 






1926.3 


4.50 


1.5 




1.0 




Fe m Sum 




13.55 


4.5 








[Sim] 

L J 


1882.7 


0.92 


0.3 


350 






Sim] 


1892.0 


11.46 


3.7 


830 








1892.0 


19.86 


6.5 


3960 






[Siml+Sim] 

L J 1 J 




32.24 


10.5 








Cm] 


1908.7 


9.5 


3.1 


875 








1908.7 


40.77 


13.4 


5500 






C m] Sum 




50.26 


16.5 






Model B a 
















Fe III UV34 


1895.5 


3.76 


1.2 




0.375 






1914.1 


10.04 


3.3 




1.0 






1926.3 


4.24 


1.4 




0.425 




Fe III Sum 




18.04 


5.9 








[Sim] 


1882.7 


1.02 


0.3 


335 






Sim] 


1892.0 


11.80 


3.9 


800 








1892.0 


17.45 


5.7 


3950 






[Sim]+Sim] 




30.27 


9.9 








Cm] 


1908.7 


9.24 


3.0 


880 








1908.7 


42.61 


14.0 


5500 






C m] Sum 




51.85 


17.0 






Model C 
















Fe m UV34 


1895.5 
1914.1 


6.33 
8.98 


2.1 
3.0 




0.7 
1.0 






1926.3 


4.50 


1.5 




0.5 




Fe III Sum 




19.81 


6.5 








[Sim] 


1882.7 


1.02 


0.3 


335 






Sim] 


1892.0 


10.12 


3.3 


800 








1892.0 


17.49 


5.7 


3950 






[Sim]+Sim] 




28.63 


9.3 








Cm] 


1908.7 


9.24 


3.0 


880 








1908.7 


4.26 


14.0 


5500 






Cm] Sum 




51.87 


17.0 







a These model measurements are also listed in Tables 3 and 4. 



